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Abstract. This paper is concerned with the three-dimensional equations of a simplified hy- 
' drodynamic flow modeling the motion of compressible, nematic liquid crystal materials. The 

\^ • authors establish the global existence of classical solution to the Cauchy problem with smooth 

0> I initial data which are of small energy but possibly large oscillations with constant state as far- 

' field condition which could be either vacuum or non-vacuum. The initial density is allowed to 

■ vanish and the spatial measure of the set of vacuum can be arbitrarily large, in particular, the 

! initial density can even have compact support. As a byproduct, the large-time behavior of the 

solution is also studied. 

Keywords. Compressible nematic liquid crystal flow; Vacuum; Large oscillations; Global well- 
^ ' posedness; Large-time behavior 

1 Introduction 

Nematic liquid crystals are aggregates of molecules which possess same orientational order 
and are made of elongated, rod-like molecules. The continuum theory of liquid crystals was 
established by Ericksen [7j and Leslie |25j during the period of 1958 through 1968, see also the 
book by de Gennes [ilj. Since then, there have been some important research developments in 
liquid crystals from both theoretical and applied aspects. When the fluid containing nematic 
liquid crystal materials is at rest, we have the well-known Ossen-Frank theory for static nematic 
liquid crystals, see the pioneering work by Hardt-Lin- Kinder lehrer [13] on the analysis of energy 
minimal configurations of nematic liquid crystals. In general, the motion of fiuid always takes 
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place. The so-called Ericksen-Leslie system is a macroscopic continuum description of the time 
evolution of the materials under the influence of both the flow velocity field u and the macroscopic 
description of the microscopic orientation configurations d of rod-like liquid crystals. 

When the fluid is an incompressible viscous fluid, Lin [27] first derived a simplified Ericksen- 
Leslie equations modeling the liquid crystal flows in 1989. Later, Lin and Liu [28^129] made 
some important analytic studies, such as the existence of weak/strong solutions and the partial 
regularity of suitable solutions of the simplified Ericksen-Leslie system, under the assumption 
that the liquid crystal director field is of either varying length by Leslie's terminology or variable 
degree of orientation by Ericksen's terminology. Recently, Wang |38j obtained the global well- 
posedness of the hydrodynamic fiow of nematic liquid crystals in the entire space under some 
small conditions. However, when the fiuid is compressible, the Ericksen-Leslie system becomes 
more complicate and there are very few analytic works available yet. The readers can refer to 
the recent works due to Morro |35] and Zakharov-Vakulenko [39] on the modeling and numerical 
studies, respectively. 

In this paper, we consider a simplified version of Ericksen-Leslie equations modeling the 
hydrodynamic flow of compressible, nematic liquid crystals in the whole spatial domain R^: 

pt + dw{pu) = 0, (LI) 
{pu)t + div(pu (g)u) + VP{p) = pAu + {p + A)Vdivn - Vd ■ Ad, (1.2) 
dt + u-Vd= Ad+\Vd\'^d, (1.3) 

where p : x [0, oo) — > 1R+ is the density of the fiuid, u : M'^ x [0, oo) — )- is the velocity field 
of the fiuid, P : x [0, oo) R+ is the pressure of the fiuid, and d : x [0, oo) (the 
unit sphere in R'^) is the unit- vector field that represents the macroscopic/continuum molecular 
orientations. The viscosity coefficients p, A satisfy the physical conditions: 

p>0, X + ^p>0. (1.4) 

The pressure P{p) is usually determined through the equation of states. Here, we focus our 
interest on the case of isentropic fiows and assume that 

P{p)=Ap^ with A>0,7>1. (1.5) 

Though the system (ll.ip -- (ll.3p is a simplified version, it still retains most of the inter- 
esting mathematical properties (without destroying the basic nonlinear structure) of the origin 
Ericksen-Leslie model for the hydrodynamics of nematic liquid crystals (cf. [6|l8|ll2p25p27j ). From 
the viewpoint of partial differential equations, ()l.iP " (|1.3p is a coupled hyperbolic-parabolic sys- 
tem with strong nonlinearities, and thus, its mathematical analysis is full of challenge, especially 
in the case that vacuum states are allowed (i.e. the density p may vanish). It is worth pointing 
out that when d = Constant, (jl.ip - (11.3p reduces to the famous Navier-Stokes equations which 
describe the three-dimensional motion of compressible viscous isentropic fiows and have been 
studied by many authors, see, for example, [9l|T7l[32l[34] , and the references cited therein. 

We shall look for the solutions, {p{x, t),u{x, t), d{x, t)), to the Cauchy problem of (jl.ip ~ (|1.5p 
with the far-field behavior 

{p,u,d){x,t) ^ {p,0,l) as |a;| — )• oo, t > 0, (1-6) 

and the initial data 

{p,u,d){x,0) = {po,uo,do){x), xeM.^, (1.7) 
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where /5 > is a given nonnegative constant, " 1 " is a given unit vector and \do\ = 1- 

The local-in-time strong solutions to the initial value or initial-boundary value problem of 
(ll.ip -( fL3]l with nonnegative initial density were studied by Ding-Lin- Wang- Wen ^ and Huang- 
Wang [20] in 1-D and 3-D spatial domain, respectively. Recently, based on the arguments 
in |18lll9j for the compressible Navier-Stokes equations, some interesting results on the blow- 
up criterion of strong solutions of (jl.ip ~ (jl.3p were obtained (see [2Tl[22l[33]). However, to the 
author's knowledge, the global existence of strong/classical solutions with vacuum is still open 
even that the initial data are suitably small in some sense. 

Recently, Huang-Li-Xin [17\ established the global existence of classical solution with large 
oscillations and vacuum to the compressible isentropic Navier-Stokes equations. Motivated by 
[17] . in this paper, we will study the global well-posedness of strong/classical solutions to the 
Cauchy problem (ll.ip -( fT77|l when the initial data are sufficiently smooth and are suitably small 
in some energy-norm. 

Before formulating our main results, we first explain the notations and conventions used 
throughout this paper. For simplicity, we set 



For 1 < r < oo and /3 > 0, the standard homogeneous and inhomogeneous Sobolev spaces are 
simply denoted as follows: 



fdx 




fdx. 



< 



l^fc.r ^ 2.'- n Z)'='^ H^ = W'''^, D'' = D^'^, = {u £ \ \\Vu\\l2 < oo 



} 



H^ = {f:R^^R\ WfWl^ = J mliOm < oo} , 



where / denotes the Fourier transform of /. 
For given initial data {po,UQ,do), we define 




where G(-) is the potential energy density given by 




It is clear that 




if 



if 



p = 



P>0, < p < 2p, 



for positive constants ci{p, p),C2{p, p) depending on p and p. 

Our main results in this paper now can be stated as follows. 
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Theorem 1.1 For given numbers Mi, M2 > (not necessarily small), p > p + 1, /3 G (1/2, 1], 
and q G (3,6), assume that the initial data (pQ,UQ,do) satisfy 



'G(po) + Po\uo\'^ + iVdop £L\ < po{x) < p, 
{po-p,P{po)-P{p))eH^nW^'i, 

uoeH^nD^nD'^, Vdo e nD^ nD^, 

Uuo\\h^<Mi, II Vdo 11^/5 <M2, \do\ = l, 
and that the compatibility condition 

1 /2 

— pAuo - (A + /x)Vdiviio + VP(po) + Vdo • Ado = Po 9 



(li 



(1.9) 



holds for some g . Then there exists a positive constant e > 0, depending only on p, X, A, 
7, p, p, Ml and M2, such that if 

Co<e, 



the Cauchy problem \1.1\ )~ [T7?^ has a unique global classical solution (p, u, d) satisfying 

0<p{x,t)<2p for all x G M^ t > 0, 

and 

'{p-p, Pip) - Pip)) G C([0, T];H^n W^'"), 
u G C([0, T]; n n L~(r, T; n L>3,9)^ 
ut G L~(r, T; n Z)^) p //^(t, T; D^), 
VdG Ci[0,T]; H"^) n L°°iT,T; H^), 
VdteCi[0,T];L^)nH^iT,T;L^) 
for any < r < T < 00. Moreover, the following large-time behavior: 



(1.10) 
(1.11) 



(1.12) 



lim (\\pi-,t)-p\\Lp+ [ p^/^\u\^dx + \\Vui-,t)\\Lr + \\Vdi-,t)\\wi,k 
holds for r G [2,6), k £ (2,6), and 



0, 



(1.13) 




if P = 0, 
if p>0. 



(1.14) 



Remark 1.1. The solution obtained in Theorem 11.11 becomes a classical one away from the 
initial time. Moreover, the oscillations of (p, u, Vd) could be arbitrarily large and the interior 
and far field vacuum states are allowed. 

Remark 1.2. When d = Constant, the system ()l.ip - ()1.3p reduces to the well-known Navier- 
Stokes equations for compressible isentropic flows. So, Theorem 11.11 improves the result due to 
Huang-Li- Xin [T^, since the compatibility condition (|1.9p imposed on the initial data is much 
weaker than the one in |17] (see also [3]). Indeed, to prove the existence of a classical solution 
of the compressible Navier-Stokes equations, the authors in ^[17] had to assume that 



— pAuo — ip + A)Vdivno + VP(/9o) = pg with g G D^, ^ -^^ 



(1.15) 
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Remark 1.3. In |10p23j. Fujita-Kato and Kato proved that the incompressible Navier-Stokes 
system is globally well-posed for small initial data in the homogeneous Sobolev spaces or 
in . In our case, since the initial energy is small, we need the boundedness assumptions on the 
-norm on the initial velocity which is analogous to the one for the compressible Navier-Stokes 
equations [E]. Note that ^ L^/{3-2/3) 6/(3-2/3) > 3 for /3 > 1/2. Thus, compared with 
the results in [101123] . the conditions on the initial velocity may be optimal under the smallness 
conditions on the initial energy. 

Remark 1.4. Recently, Wang [37[|38j proved that the heat flow of harmonic maps (11. 3p (with 
= 0) and the incompressible liquid crystal flow are globally well-posed provided that ||V(io||L3 
and IIuoIIls + llVdollLS are sufficiently small, respectively. In view of these results in [371138]. 
the conditions imposed on the initial director field do may also be optimal under the smallness 
conditions on the initial energy. 

We now comment on the analysis of this paper. For large initial data satisfying (jl.Sp 
and (jl.9p . one can utilize the Galerkhin approximation method to construct the local classical 
solutions in a similar manner as that in [21] (see, Proposition 15.11 below). So, to extend the 
classical solution globally in time, we need some global a priori estimates on the solutions 
(/3, u, d) in suitable higher norms. To do so, we notice that (ll.ip - (ll.3p are indeed a coupled 
system between the Navier-Stokes equations and the equations for heat flow of harmonic maps, 
so that, we shall make use of some ideas developed in [HKlT]. However, compared with the 
compressible Navier-Stokes equations, some new difficulties arise due to the additional presence 
of the liquid crystal director field d and the weaker compatibility condition ()1.9I) (cf. (jl.lSp ). 
Especially, the super critical nonlinearity |V(ipd in the transported heat flow of harmonic map 
equation (|1.3p and the strong coupling nonlinear term Vd • Ad in the momentum equations (jl.2p 
will cause serious difficulties in the proofs of the time-independent global energy estimates. 

As that in [T7], it turns out that the key issue in this paper is to prove both the time- 
independent upper bound for the density and the time-dependent higher norm estimates of the 
solutions (p, M, d). For this purpose, as usual we start with the basic energy estimate (see Lemma 
13. ip . To overcome the difficulties induced by the nonlinearities | Vdpd and Vd • Ad, we succeed 
in deriving an estimate on the spatial L^-norm of the gradient of d and the estimates on the 
second-order spatial derivatives of d, which are indeed suitably controlled by the initial energy 
Co and the ij^-norm of Vd (see Lemmas 13.21 13. 3p . Then, basing on these estimates, we can 
utilize the techniques in [141I151 [T7] to obtain an estimate on the spatial L^-norm of the velocity, 
and to carry out some careful initial-layer analysis which is concerned with the elegant estimates 
on the gradient of the velocity, the second-order and third-order derivatives of the director field, 
and the material derivatives of the velocity as well (see Lemmas l3.4H3.7p . It is worth pointing 
out that similar to that for the compressible Navier-Stokes equations, the effective viscous flux 
and the vorticity (see p.3p for the definition) play a very mathematically important role in the 
entire analysis, which are instrumental in controlling || Vii||_LP (2 < p < 6) by the L^-norm of the 
gradient and the material derivative of the velocity (see Lemma 12.21) . With these estimates, we 
then can obtain the desired estimates on both -L^(0, min{l, T}; L°°) and L^/'^(min{l, T}, T; L°°) 
norms of the effective viscous flux, and thus, it follows from Zlotnik's inequality (see Lemma 
12. 3p that the density admits a uniform-in-time upper bound which is the cornerstone for the 
global classical estimates of the solutions. 

The next main step is to estimate the derivatives of the solutions. Indeed, to achieve these, 
we first apply the Beale-Kato-Majda type inequality (see Lemma 12. 4p to prove the important 
estimates on the gradients of the density and velocity by solving a logarithm Gronwall inequality 
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in a similar manner as that in |17tll9j. As a result, one can easily obtain the L^-estimates for 
the second-order derivatives of density, pressure and velocity. However, due to the weaker 
compatibility condition in (II. 9p (cf. (I1.15P ). the method used in [TTllll] cannot be applied any 
more to obtain further estimates. Motivated by |161I26|. instead of the L^-method, we succeed 
in achieving these classical estimates by proving some desired L'^-estimates (3 < g < 6) on 
the higher-order time-spatial derivatives of the density and velocity, basing on some careful 
initial-layer analysis (see Lemmas I4.4t44.6| l . 

The rest of this paper is organized as follows. In Sect. 2, we first collect some known 
inequalities and facts which will be frequently used later. In Sect. 3, we derive the time- 
independent (weighted) energy estimates of the solutions and the key pointwise upper bound of 
the density, which will be used to study the large-time behavior. In Sect. 4, we establish the 
time-dependent estimates on the higher-order norms of the solutions, which are needed for the 
existence of classical solutions. Finally, the proof of the main result (i.e. Theorem 1.1) will be 
done in Sect. 5 



2 Auxiliary lemmas 

In this section, we recall some elementary inequalities and known results which will be used 
frequently later. We start with the well-known Sobolev inequalities (see, for example, |24j). 

Lemma 2.1 For 2 < p < 6, l<g<oo and 3 < r < oo, there exists a generic constant C > 0, 
depending only on q and r, such that for f € and g ^ L'^ r\ D^'^ , we have 

II/IIl. < c||/||g-^)/(^^)||v/||(.r^)/(^^\ (2.1) 

Ml^ < cii^iiir'^/^'^+'^^'^-'^^iiVffiii':/^''^^''^''-')^ (2.2) 

As that for compressible Navier-Stokes equations (see, for example, [I11[I71[32]), the con- 
nections among the effective viscous flux, the vorticity and the physical quantities will play an 
important role in the entire analysis of the present paper. So, to be continued, we set 

F = (2/i + A)divn- (P(p) -P(p)), uj = Vxu, M{d) = VdQVd - ^\Vd\%, (2.3) 

where F is the so-called effective viscous flux, uj is the vorticity, I3 is the 3x3 unit matrix, and 
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VdQVd = Y, did^djd^ 



k=l 

Then, it follows from Lemma |2. II and the standard L^-estimates of elliptic system that 



Lemma 2.2 Let {p,u,d) be a smooth solution of \1. 1] )- (T77^ on x (0, T]. Then there exists 
a generic constant C > 0, which may depend on fi and A, such that for any p G [2, 6], 

\\VF\\lp + WVloWlp < C {WpuWlp + W\^d\\V^d\WLp) , (2.4) 

WHl^ + W^Wl^ < C (WpuWl^ + WNd\\V^d\WL2) , (2.5) 

and 

WVuWl^ < C {WpuWl^ + II^'(P) - P{P)Wl<^ + IllVdllV^dlll^.) . (2.6) 
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Proof. Indeed, due to (ll.2p . one has 

pu = VF - nV X u - div(M((i)), 

and hence, 

AF = dw{pu) + divdiv(M((i)), /xAw = V x (pu + div(M((i))), (2.7) 

where f — ft + u ■ V/ denotes the material derivative. Thus, an apphcation of the standard 
L^-estimate of elhptic system leads to (j2.4p . which, together with (|2.ip . gives (|2.5p . Using (j2.3p . 
(|2.5p and the standard L^'-estimate, we obtain (|2.6p . □ 

To prove the uniform-in-time upper bound of density, we need the following Zlotnik in- 
equality, whose proof can be found in |40j . 

Lemma 2.3 Assume that the function y £ W^'^{0,T) solves the ODE system: 

y' = g{y)+b'{t) on %T], y(0) = yo, 

where b € W^^'^{0,T) and g G C(R). // 5(00) = -00 and 

bit2) - Kti) < No + N,{t2 - ti) (2.8) 

for all < ti < t2 < T with some Nq > and Ni > 0, then one has 

y{t) < max{yo,r} + iVo < +00 on [0,T], (2.9) 

where ^* G M is a constant such that 

9{0<-Ni for e>r- (2.10) 

Finally, we recall the following Beale-Kato-Majda type inequality (cf. [Il[l9]), which is an 
essential tool for the estimates of the gradients of {p,u). 

Lemma 2.4 For q £ (3, 00), there exists a constant C{q) > such that for all Vu G n D^''^, 
||Vu||loo < C (lldivnllioo + ||V x n||L°o) In (e + \\V'^u\\li) + C||Vn||i2 + C. (2.11) 

3 Time-independent lower-order estimates 

This section is devoted to proving the time-independent (weighted) energy estimates and 
the uniform upper bound of density. Assume that (p, u, d) is a smooth solution of (ll.l[ )- (jl.7p on 
X (0, T) with some positive time T. To estimate the solutions, we set 

a(t)^{l,t} 

and define the following functionals: 

Ai{T)^ sup a{\\Vu\\l, + \\V^d\\l,) 

0<t<T 

+ ^"^^1 (||p'/'n|li2 + \\Vdt\\l2 + llV^dlll^) dt, (3.1) 
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A2iT)^ sup am\p'/\\\l, + \\Vdt\\l2 + \\V''d\\l2 

0<t<T 



+ / am\Vu\\l, + \\du\\l2 + \\V^dt\\l2)dt, 
AsiT)^ sup \\Vd\\l3+ r [ \Vd\\V'^d\'^dxdt 

0<t<T Jo J 

A,{T)^ sup CT(3-W4(||v^||2^ + ||v2d||2,) 
0<t<T 



(3.2) 
(3.3) 



+ 



a 



(3-2/3)/4 



p^'^u\\l2 + \\V^d\^i2)dt 



and 



A(T) ^ sup / p\ 

0<t<T J 



u\ dx, 



(3.4) 
(3.5) 



where /3 E (1/2, 1] and v = 't;t + n • is the material derivative. 

We shah prove the fohowing key a priori estimates on the solutions {p,u,d). 

Proposition 3.1 For given numbers Mi,M2 > 0, p > p + 1 and f3 £ (1/2,1], assume that 

'G{po) + po\uo\'^ + \Vdo\^ G L\ 

< < infpo < sup/9o < p, \do\ = 1, (3-6) 
Muo\\hp<Mi, ||Vdo|lij/3 < Ms. 
Suppose that {p, u, d) is a smooth solution of Ill.l\) - p77\ ) on M'^ x (0, T] satisfying 



0< inf p{x,t) < sup p{x,t)<2p, 

(x,t)eR3x[0,T] (x,t)eK3x[0,T] 



^i(r) + A2(T) < 2Co'/', AsiT) < 2Co'«, A^{a{T)) + A^iaiT)) < 2Co'°, 



So 



where 



A 2/3-1 
9P 



G (0,1/9] for /3g (1/2,1]. 



(3.7) 



(3.8) 



Then there exists a positive constant e > 0, which depends only on p, X, Ajj, p, p, 13, Mi, and 
M2, such that 



provided 



< inf p{x,t) < sup p{x,t) < jp, 

{x,t)GlR3x[0,T] (a;,t)GR3x[0,T] 

AiiT) + A2iT) < Cy", A^{T) < Co^°, A^iaiT)) + A^iaiT)) < C',\ 



Co<e. 



(3.9) 



(3.10) 



The proof of Proposition 13.11 is based on a series of lemmas and is postponed to the end 
of this section. For notational convenience, throughout this section we denote by C or Ci 
(i = 1, 2, . . .) the generic positive constants which may depend on p. A, A, 7, p, p, j3. Mi, and M2, 
but are independent of T. We also sometimes write C(a) to emphasize the dependence on a. 

We begin with the following estimates. 
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Lemma 3.1 Let {p,u,d) be a smooth solution of / fj. j|) - |777| j on x (0,T] satisfying |g. 7\ ). 
Then there exists a constant ei > 0, depending only on n, X, A,^, p, p, Mi, and M2, such that 



and moreover, 



provided Cq < ei- 



sup f ip\u\^ + G{p) + \Vd\^)dx 

<t<T J 

+ r + WdtWh + l|V'd|li2) dt < CCo, (3.11) 
Jo 

r (ll^llle + llVnlli^ + llV^dll^,) dt < CCf\ (3.12) 
Jo 



Proof. Multiplying (|1.1|) and (jl.2p by G'{p) and u in respectively, integrating by parts and 
adding them together, by (j2.ip we have 

^ y" (^G(p) + ^H^^l')^^^ + j (/i|Vn|2 + (^ + A)(divn)2)dx 
= - j{Vd- Ad) ■ udx < C\\u\\l6 ||Vd||i3 II V2(i||i,2 

<fl|Vn||i.+C7||Vd||i3||V^d||i.. (3.13) 
Using the fact that |(i| = 1 and integrating by parts, we infer from (jl.3p and (j2.ip that 

ivd|2dx+ y {\dt\'^ +\v^d\^)dx 



d 

di 



j \dt-Ad\'^dx = j \u-Vd-\Vd\ 



^dl^ dx 



< 



c(IHIi.||vd||i3 + ||vd||i4) 



<C||Vd||i3(||V^z||i. + ||v2d||i.), 
which, together with ()3.13p and (j3.7p . gives 

I / [G{p) + lp\u\' + \Vd\'yx+(^^\\Vu\\h + \\d,\\^^^ 
< C\\Vd\\ls {\\Vu\\l, + llV^dlli.) < CiCo^'"/^ (llVnIli. + . (3.14) 

Thus if Co is chosen to be such that 



Co < ffi = min < 1, 



^ n3/(25o) / J \3/(25o)- 



4CJ 'UcJ 



then integrating (I3.14p in t over [0,T] immediately leads to (13. lip . 
To prove (l3T2]l . we utihze ([21]), dST]) and (l3TT]) to deduce that 



^(IHli6 + ||Vn||i. + ||V24t,)di 
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< C / (||Vn||i2 + llV^dll^,) dt+ G (||V^/||^2 + llV^dlli^) dt 

Jo J'^iT) 

<C sup U'-'^ym\Vu\\l, + \\V''d\\l,)) / 

0<t<a{T) ^ ^ Jo 

+C sup a{\\Vu\\l2 + \\V'^d\\l2) r {\\'^u\\l2 + \\V^d\\l2)dt 

a(T)<t<T Ja(T) 



since /3 S (1/2, 1] and 5o G (0, 1/9]. The proof of Lemma [3TT] is therefore completed. □ 

The next lemma is concerned with the estimate of A^{T). 

Lemma 3.2 Let {p,u,d) be a smooth solution of lll.l\) - [T7^ on x (0,T] satisfying |3. 7\ ). 
Then there exists a constant £2 > 0, depending only on fi, X, A,j, p, p, f3, Mi and M2, such that 

A3{T)+ r \\Vd\\%dt<c'Q°, (3.15) 
Jo 

provided Cq < £2- 

Proof. Operating V to both sides of (jl.3p gives 

Vdt - VAd = V{\Vd\'^d) - V(n • Vd). (3.16) 
Thus, multiplying (j3.16p by 3\Vd\Vd and integrating by parts over M'^, we obtain 
d 



dt' 



;||Vd||i3+3 J (\Vd\\V^d\^ + \Vd\\V{\Vd\)\^^dx 
<cj {\Vd\^\V^d\ + \u\\Vd\'^\V^d\) dx 

<\j \Vd\\V^d\''dx + C[\\Vd\\l, + \\Vu\\l4Vd\\l,^,). (3.17) 
To deal with the right-hand side of (|3.17p . we first use Lemma l2.1l to get that 

\\vdu. < c||vd||^/^||vd||J/^ iivdiL^/. < c\\vd\\'ff\\vd\\'ff, 
livdili, = |||vd|3/2||2, < c||v(|vd|3/2)||2^ < c\\\vd\y^vM\h, 

so that, by virtue of (|3.7p and Cauchy-Schwarz inequality we find 

R.H.S. of (3.17) < (1 + C\\Vd\\l^) |||Vd|^/2y2^||2^ ^ ^||^^||3^||^^||4^ 

< (1 + CiCf'''^) |||Vd|i/2v2d||2, + CCl^WVufL,. (3.19) 



Thus if Co is chosen to be such that 

Co<£2,i = min{ei,C7r'/(''°)}, 
then we obtain after putting ()3.19p into (j3.17p and integrating it over (0, T) that 

sup ||Vd||i3+ r (\\Vd\\y + \\\Vd\'''V^d\\l,)dt 

0<t<T Jo ^ ' 
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(3.18) 



J \\Vu\\l2dt 

< CCo^(2/3-l)/(4/3) ^ ^^35o < (^2Co^'^^ (3.20) 



< C\\Vdo\\ls+C&o" 



where we have used ([3J2]l . (fXTH]) . (|2J]l . and the Sobolev embedding inequahty ^ L^/i^-'^P). 
Therefore, by choosing Co sufficiently small to be such that 



Co<e2 = min{e2,i,^^2~'^°} 



we immediately obtain the desired estimate in (j3.15p from (|3.20p . □ 
The following initial-layer estimate of d is crucial for the subsequent analysis. 

Lemma 3.3 Let [p, u, d) he a smooth solution of ( fi.i| )-( f777| j on M'^ x (0, T] satisfying {3. 7| j. 
Then there exists a positive constant e^,, depending only on p, X, A,^, p, p, P, Mi and M2, such 
that for any £ [0, 1], 

sup U-'WV'dWl,) + r a^-' {\\Vdt\\l2 + W^mi,) dt < CIlVdoIlL, (3.21) 

<t<T ^ ^ Jo 



0<t<T ^ ^ Jo 

provided Cq < £3 . In particular, if Cq < £3 , then 

{o\w'd\\i^ + r 

0<t<T 



sup (cxIlV^dll^^) + r a{\\Vdt\\l2 + \\V^d\\l2)dt<C\\Vdo\\l2 <CCo. (3.22) 

<t<T Jo 



Proof. Let {p, u, d) be the smooth solution of ()l.iP " (|1.7p on R"^ x (0, T]. Consider the following 
Cauchy problem of linear parabolic equations: 



vt - Av = -u-Vv + dVd -.Vv, v{x,0) = vo{x), (3.23) 



where A : B = Y^fj^i aijbij for A = (aij)3x3 and B = (bij )3x3- 

By (j2.ip and Holder inequality, we easily deduce from (|3.23p that 

j^W\i2 + {\\vt\\h + \\vMh) 

= J \vt — Av\'^dx = J \u ■ Vv — dVd : Vvl^dx 

< C\\Vu\\12\\Vv\\l2\\Vv\\lc> + C||Vd||2g||Vt;||^//||V7;||^i^ 

< IW'^Mh + C {\\Vn\\l2 + \\Vd\\%) \\Vv\\l2, 
which, together with Gronwall's inequality, (j3.12p and (|3.15p . gives 

sup \\Vv\\l2+ [ {\\vt\\l2 + \\VMh)dt<C\\Vvo\\l2. (3.24) 

0<t<T Jo 

Applying V to both sides of (|3.23p and taking the L^-inner product, we get 
|l|V^-lli. + (l|V^;.||i. + ||V=^Hli.) 
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< J {\V{dVd:Vv)\'^ + \V{u-Vv)\'^)dx = I. (3.25) 
Using Lemma |2. 11 (|3.15p and Cauchy-Schwarz inequality, we know that 

/ < [(iVdp + |n|2) iV^wl^ + (iV^dp + |Vd|^ + |Vu|2) jVul^] dx 

< C\\Vd\\l,\\V'^v\\l, + C||u||2,||V2|;||^2||V27;||i6 
+C\\Vv\\U iW^mh + llVdlliallVdllie + llVulli.) 

< Q + CaCo''"/^') ||v3t;||2, + C\\V\\\l, {\\V'd\\l, + \\Vu\\l,) . (3.26) 



Thus if Co is chosen to be such that 

-'2 



Co < £3 = (4C2^~=^/(25o) 



then substitution of (I3.26P into (I3.25P results in 

^llV^^lli. + (||V^;,||i. + llV^^lli.) < C\\VMh iW^'dWh + llV^lli.) , (3.27) 
which, together with Gronwall's inequality and p.l2p . implies that 

sup IIV^t;!!^^ + / {\\Vvt\\l2 + \\V^v\\l2)dt<C\\V\o\\l2. (3.28) 

0<t<T Jo 



On the other hand, multiplying (j3.27p by a and integrating it over (0,T), we infer from 
([3:2i]l that if Co < £3, then 



sup {(7\\V\\\l2)+[ a {\\Vvt\\l2 + \\V'^v\\l2) dt 

<t<T Jo 



0<t<T 

t-T 



<C I \\V^v\\l2dt <C\\Vvo\\l2. (3.29) 





Since the solution operator vo 1— ?• v{-,t) is linear, one may apply the standard Riesz-Thorin 
interpolation argument (see 12j) to (j3.28p and (|3.29p to get that for any 9 £ [0, 1], 

sup (t'-'\\VMl2)+ r t'-' {\\Vvt\\l2 + \\VMh)dt<C\\Vvo\\%e, (3.30) 

0<t<T ^ ^ Jo 

with a uniform constant C independent of 9. Thus, choosing vq = do so that v = d, then (I3.2ip 
follows from (|3.30p directly. □ 

By Lemmas I3.1H3.31 we can now derive preliminary bounds for Ai{T) and A2{T). 

Lemma 3.4 Let {p, u, d) he a smooth solution of / fi.i|) -( (777| ) on E'^ x (0, T] satisfying {3. 7\ ). 
Then there exists a constant £4 > 0, depending only on p, X, A,j, p, p, f3, Mi and M2, such that 

Ai{T)<CCo + C r a''{\\Vu\\\, + \\P{p)-P{p)\\\,)dt, (3.31) 
Jo 

and 

A2(T) < CCq/^^^^" + CAi{T) + C / a'^\\Vu\\lidt. (3.32) 

Jo 

provided Co < £4. 
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Proof. Multiplying (II. 2|) by aut and integrating by parts over give 

= ^ y fj ((P(p) - P(p)) divn + M{d) : Vu) dx 

+a J pu ■ Vn • iidx ~ ^ J i-^ip) ~ P{p))t divndx ~ J -^(^)t • Vndx 

+a' J Q (/u| Vn|2 + + A)(divn)2) - (P(p) - P{p)) divu - M(d) : Vu^ dx 

"i^o + E^- (3.33) 

1=1 

We are now in a position of estimating the terms on the right-hand side of ()3.33p . First, by 
Holder and Cauchy-Schwarz inequalities we easily see that 

< la\\p^/\\\l2 + Ca^\\Vu\\li + C\\p^/^u\\l4. (3.34) 
8 

To deal with I2, we first deduce from (jl.ip that 

{P{p) - P{p)\ + n • V (P(p) - P{p)) + 7P(p)divu = 0, (3.35) 
which, together with the effective viscous flux F in ()2.3p . yields that 

^2 = crJ 'yP{p)idivufdx + a j u ■ V {P{p) - P{p)) divudx 

< C\\Vu\\l2 + Ca^\\P{p) - P{p)\\U + C<y j \P{p) - P{p)\ (|Vn||F| + \u\\V F\) dx 

< C\\Vu\\l, + Ca^Pip) - P{p)\\l, + Ca\\Pip) - P{p)\\l4'^Al4'^ 

< C\\Vu\\l, + Ca^Pip) - P{p)\\h + Ca\\Vu\\L2 (Wp'^^iih^ + \\\V d\\V^ d\\\ ^2 



< IcrWp'^'^h + C\\Vu\\l, + Ca\\Vmi, + Ca^\\P{p) - P{p)\\h, (3.36) 

where we have also used (j2.ip . (|3.1ip . (|3.15p and Cauchy-Schwarz inequality. 
It follows from (j2.ip . (j3.15p and Cauchy-Schwarz inequality that 

/3 < Ca\\Vd\\]^,^\\Vd\\]^,^\\Vdt\\L2\\Vu\\L, 

< C\\V^d\\l2 + Ca\\Vdt\\l2 + Ca^\\Vu\\l4, (3.37) 

and finally, it is easily seen that 

h < ca' imp) - pmh + \Mh + iivdiii4) 

< Ca' (Co + llVnlli^ + C^^'^^WV^dWl^') . (3.38) 
13 



Due to (IXTTD . (fXT5]) and ([322]), we have 

h < ^cxIlV^lli, + Ca {\\P{p) - P{~p)\\l, + ||Vd||i3 \\V^d\\l,) 

<^(y\\yu\\l2+CCo, (3.39) 

so that, putting (|3.34p and ()3.36p - ()3.39p into ()3.33p and integrating it over (0,r), we obtain 

rT 

sup (frllVnllla) + / aWp^/^iiWl^dt 
0<t<T Jo 

rT 

|2 I |iv72j||2 , _|iv73j||2 , _||V7J l|2 



<CCo + C {\\Vu\\l, + \\V^d\\l, + a\\V'd\\l, + a\\Vdt\\l,) dt 
Jo 

+C ra'{\\Vn\\l + \\P{p)-P{p)\\Udt + C ^ Wp^MUt 
Jo Jo 

<CCo + C [ am\Vu\\l, + \\P{p)-P{p)\\h)dt + C [ \\p'/^u\\l,dt, (3.40) 
Jo Jo 

provided Co < £3. Here we have also used ()3.1ip and (|3.22p . 

To estimate the last term on the right-hand side of (|3.40p . we observe from (j2.ip . (|3.7p and 
(I31T]1 that 

' \\p'/Mh < CWp^MhlNuWl, < C\\Vu\\l, for 0<t<a{T), 
Wp^^Mh < CWp^Ml^ \\Vu\\l, < C\\Vu\\l, for a{T) <t<T, 
which, combined with (|3.11|) . gives 



[ \\p^/^u\\l4dt <C [ \\Vu\\l2dt < CCo. 
Jo Jo 



This, together with (lOOll and finishes the proof of ([MI])- 

To prove (I3.32p . operating a'^ii^ [dt + div(n-)] with m > to both sides of the j-ih. equation 
of ()1.2p and integrating by parts over M'^, we obtain after summing up that 

— {p + A) y (j™tf' [djdiYUt + div(n(9jdivu)] dx 
= ^(j^-Vy p\u\^dx - j a'^ii^ [djPt + div(n5jP)] dx 

4 

+ / a"'M{d)t : Viidx + j (j™divM(d) • (n • Vu)dx = ^Ri. (3.41) 

i=l 

After integrating by parts and using Cauchy-Schwarz inequahty, we easily see that 

dx 



Li^-pj a'^v? [Au^' + div(uAu-') 

= (I v„P-a.„^*(„.v„^)+%.^„'A„.).. 
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>T^<j'"l|Vil|li.-C,T"1Va||*., 



where and in what fohows, we use the Einstein convention that repeated indices denote the 
summation with respect to those indices. Similarly, 

L.-(, + A)/,"„Ma,d,v„, + d„(„a,d™)i<ix 

>(li + X}a'^\\iim\\l, - ^J<t"'||Vu|||j - C(t'"||Vu||'.. 

o 

Thus, the left-hand side of ()3.4ip can be bounded from below as follows: 

L>~ {o-^Wp^'^uWh) + f cT™||Vn||i. - Ca^\\Vu\\\,. (3.42) 
In view of (j3.35p . we have by integration by parts and (j3.7p that 



i?2 = / ct" ydkU^v^djPip) - divii (div(P(/))n) + (7 - l)P(p)divn)J dx 
(dkii^dju'^Pip) + (7 - l)P(p)(divn)(divn)) dx 
< ^'^"l|Vn||i2 + Ca"^\\Vu\\l,. (3.43) 

o 

Using (HU, ([22]) and (f3l3D . we get 

R3<Ca^\\^d\\L4^dth4'^uh^ 

< ^^'"llVnlli, + CiC^^'^'^'a^WV^dtWl^, (3.44) 



and similarly. 



i?4 < Ca'^\\VdUoo\\V^d\\]^,^\\V^d\\];,^\\u\\L4Vu\\L2 



< |a-||Vu||i. +C7a-||V2d||i.||V3d||2,||Vn||i. 

< f ^'"llVnlli. + Ca^ {\\Vn\\l, + \\V'd\\l,) WV^l,. (3.45) 
Thus, putting ([332]) -([335]) into ([3311), we know that 

< Cma''"^a'\\p^/^u\\l2 + Ca"'\\Vu\\l2 + Ca'''\\Vu\\l4 

+Ca^ (llVnIli. + \\V'd\\\2) llV^dlli.. (3.46) 

Next, we estimate ||V^(it||j;^2. To do so, noticing that 

dtt-Adt= (iVdpd-n-Vd)^, 
we obtain from direct computations that 

^ j a"'\Vdt\^dx + j a"'{\dttf + \V^dt\^)dx-ma'^^^a' J \Vdt\^dx 
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<C J a'''\u\'^\Vd\^dx + C J a'''\u\'^\Vu\^\Vdfdx + C J a"'\u\^\Vdt\'^dx 

+C j a'^\Vd\'^\Vdt\'^dx + C I a'^\dt\'^\Vd\Ux^^Ji, (3.47) 

where we have used the facts that = 1 and ut = u — u ■ Vn. 

For the first term on the right-hand side of (j3.47p . it follows from (j2.ip and (j3.15p that 

Ji < Ca'^WiiWleWVdWls < C2Co^°/V'"||Vtt||i2. (3.48) 
Using dan, ([22]), (i3TT]l and (f3l3D . we see that 

J2<Ca'''\\u\\l4Vd\\le\\Vu\\l, 

< Can\yu\\l4Vd\\le {WpiiWh + II^P) - nml^ + ll|Vd||V2d|||2,) 

< Can\yu\\l4Vmi, [\\p'/'u\\h + Co^' + llVdlliallV^dllia) 



< Ca^ {\\Vu\\l + \\Vd\\l,) [\\p'/'u\\h + +Can\Vu\\l4Vmi., (3.49) 

and using (|2.ip and Cauchy-Schwarz inequality, we have 

5 

Y.J^< Ca^ (llnllie + ||Vd||i«) \m\\ls + Can\dt\\UVd\\% 



i=3 



< Ca^ (llVnIli, + \\V^d\\l,) \\VdtU4VHt\\L^ + Ca^-\\VM\U\^ dt\\l. 



which, combining with ()3.47p - (|3.49p gives 
d 



KllVdtlli.) + {WV'dtWh + WduWh) - C2Co^*/V™||Vn||i. 
<Cma'^-'a'\\Vdt\\l2+Ca^\\Vu\\l4V^d\\l, 

+Ca^ (llV^lli. + llV^dlliO (Wp^'MI. + ||Vd,||i. + \\V'd\\l,) . (3.50) 



Thus if Co is chosen to be such that 

Co < 54 ^ min {.3, (2Ci)-3/(2^«), (2C2)-='/(2^°)} , 

then, choosing m = 2 in (j3.46p and (|3.5Up . adding them together, and integrating the resulting 
inequality over (0,T), by p.7p and (|3.1ip we obtain 

sup a^(\\p'/Mh + \\Vdt\\l, 
0<t<T ^ 



+ / <7m\vu\\i,+\\du\\i2 + \\v'dt\\h)dt 

Jo 

<C [ a^\\Vu\\l2dt + C [ a^\\Vu\\l2\\V^d\\l2dt 
Jo Jo 
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3 j||2 
L2 



+C a' (llVnIli. + llV^dllt.) + ||Vd,||i. + \\V'd\\ 

+C j\ [\\p^/^u\\l2 + llVdtllia) dt + C j\^\\Vu\\\,dt 



dt 



< CCq^"^^^^' + C7^i(r) + C7 / a2||Vu||i4dt, (3.51) 

Jo 

since it follows from ()3.7p and ()3.12p that 

(llVnIli. + llV^dll^,) (llpi/^nlli. + llVdtlli. + llV^dlli.) dt 

'u\\l. + + llV^dll^,) (llVnIlt. + dt 



< sup (T^ 
0<t<T 

a/2+2<5o 



< CC^^°A2{T) < 

Moreover, it follows from (|3.16p and the standard L^-estimate of elliptic system that 

llV^dlli^ < C (||Vdi||i2 + iivdilie + livdv^dii^^) 



+C7 (||VwVc?||i2 + \\uV'^d\\i^2 

1^2 II V ^11^2 , 



< c ( ||Vdi||^2 + iiv^diii^ + ivHflh^^dfi^ 



so that, 



llV^dlli^ < CllV^dlli^ (llV^dlli^ + llV^lli^) + C||Vd,||i2. (3.52) 
Therefore, keeping in mind that /? > 1/2, we infer from (|3.7p and p.5ip that 

sup (a^llV^dlliO 

0<t<T 

<C sup (a2||Vdt||i2)+C7 sup (^^llV^dlli^dlV^lli^ + llV^dlli^)) 

0<t<T 0<t<T 

^^i/2+25o ^ c'Ai(r) + C7 / a'^\Vu\\\^dt 

Jo 



< 



+C sup {a\\vM\U (^^'-'^^/MllV^^lli^ + llV^dlli^))' 

0<t<(T{T) 0<t<o-{T) ^ ^ 

+C sup (a||v2d||2,) sup {<y{\\Vu\\l, + \\V'd\\l,)f 

o-{T)<t<T (T{T)<t<T 

rT 



< CCl^^^^^" + CAi{T) + C [ a^\\Vu\\l,dt, 

Jo 



which, together with (|3.5ip . finishes the proof of (|3.32p . □ 

The next lemma plays an important role in the proof of the uniform upper bound of p. The 
ideas of the proof are motivated by the ones in |15yi7]. 
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Lemma 3.5 Let {p, u, d) he a smooth solution of / fj. j|) ^ |777| j on x (0, T] satisfying {3. 7\ ). 
Then there exist positive constants £5 and C depending on /x, A, ^4, 7, p, p, /3, Mi and M2, such 
that 

sup (i^~'^||Vu||^2) + / t^-^\\p^/^u\\l2dt<C{Mi,M2), (3.53) 
t<am ^ ^ Jo 



0<t<a{T) 

and 



sup t^-l^ (\\p'^/^u\\l2 + \\Vdt\^^- 
o<t<cr{T) 

ra{T) 



lL2 ^ II vafii^a^ 

+ / t^-^ {\\V^l2 + \\du\\l2 + \\V^dt\\l2)dt<C{Mi,M2), (3.54) 
provided Co < £5. 

Proof. The proof of (I3.53P is motivated by |151ll7j. For a fixed smooth solution {p,u,d), we 
define the hnear differential operator jC acting on the functions if : x [0, 00) — >• by 

{Cwy = purj- + pu • Vw^ — [p-Aw^ + {p + X)divwxj) 

= pw^ - {pAw^ + {p + X)dwwx^) , j = 1,2,3, (3.55) 

where w = wt + u ■ Vw. We thus define wi, W2 and t^a by 

iZwi = 0, wi{x,{)) = wiq{x), (3.56) 

Cw2 = -VP{p), W2{x,Q) = G, (3.57) 

and 

Cw^ = -d\w{VdQVv - -Vd:Vvlz), W3(a;,0)=0, (3.58) 

where v = v{x,t) is the solution of (|3.23p . 

The estimates of wi,W2 are similar to those in |151 ll7j. For the reader's convenience, 
we reproduce the proofs here. Multiplying (|3.56|) and (|3.57|) by wi and W2, respectively, and 
integrating them by parts over x {{),a{T)), one easily deduces from (|3.7p and (|3.1ip that 

/rcr{T) 
p\wi\^dx+ ||Vu;i||^2(it < C||u'io||i2, (3.59) 

Jo 

and 



sup / p\w2\ dx + ||Vw2||i,2di < CCq. 

0<t<a{T) J Jo 



(3.60) 



To estimate W3, multiplying ()3.58p by W3 and integrating it by parts over M'^ x (0,(t(T)), 
by (EU and (I3T5I1 we get 

sup / p\w3\ dx+ ||V'«;3||^2iit 

0<t<a{T) J Jo 
fMT) 

<C \\yd\\LA\^v\\L4Vw3\\L2dt 
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1 f^iT) ra{T) 

JO Jo 
which, combined with (j3.24p . yields 

ra{T) 

I p\w3\^dx+ I 

0<t<a{T) 



//.ct(T) 
p\w3fdx+ \\Vw3\\l2dt<C\\Vvo\\l2. (3.61) 

. . Jo 



Next we estimate ||Vt(;j||j;^2 (i = 1,2,3). Thanks to (j3.7p . an application of the standard 
L^-estimate of elliptic system to (|3.56p shows 

||Vw;i||i6 < C\\V^wi\\l2 < C\\p^/^wi\\l2. 

So, multiplying (j3.56p by wu in and integrating by parts, we infer from (13. 7p that 

j {fi\Vwi\'^ + + X){divwif) dx + j p\wi\'^dx 

pwi ■ {u ■ V'Wi)dx 



1 d 
2di 



1/2 / ^ 1/3 



<C [ I p\wi\^dxj [ I p\u\'^dx] llVwill^e 
1/3 



<c(^Jp\ufdx^ J p\wifdx 
< CiCo°/^ / p\wi\'^dx (3.62) 



for any t G (0,(t(T)). Thus if Cq is chosen to be such that 

Co<e5,i=niin{£4,(2Ci)-3/^°}, 
then it follows from ()3.62p that 

\\Vwi\\i2 + I 

0<t<cr{T) 



sup ||Vu;i||^2+/ \\p'^^^wi\\l2dt<C\\Vwio\\l2. (3.63) 
:t<o-(T) Jo 



On the other hand, multiplying (I3.62p by t and integrating it over (0,(t(T)), by (I3.59P we 
see that if Co < £5 i, then 



f f 

sup (i||Vt(;i 11^2) + / t p\wi\ dxdt 

t<a{T) Jo J 

<C ||Vw;i||i2dt<C||t(;io||i2. (3.64) 

Jo 



Since the solution operator wiq i— )• 'Wi{-,t) is linear, by the standard Riesz-Thorin interpo- 
lation argument (see j2]), we conclude from (I3.63P and (I3.64P that for /3 G (1/2, 1], 

sup t^-^||Vu;i||i2 + / t^~^\\p^^^Mhdt<C\\ww\\l^. (3.65) 
t<crm ^ ^ Jo 



0<t<a{T) 
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In order to estimate ||Vtt;2||L2, analogous to the proof of Lemma [221 we set 

F ^ (2/i + X)diYW2 - {P{p) - Pip)). 
Then it follows from (|3.57|) that 

IIVFIU2 < C\\pw2\\l2, \\Vw2\\l^ <C{\\pw2\\l^ + \\P{p) - P{p)\\l<^) ■ (3.66) 
Now, multiplying (|3.59p by W2t in and integrating by parts, by (|3.35p we obtain 

■ j {p\Vw2\'^ + {p + X){dww2)'^)dx + j p\w2\^dx 

j {P{p) - P{p)) dww2dx - J {P{p) - P{p)) u ■ Vdivw2dx 



1 d 
2di 



d 
It 



+ j {{-i -l)P{p)+P{p)){d\wu){di\W2)dx + j pw2 ■ {u ■ V W2) dx 



~ di 



Io + ^I^. (3.67) 



i=l 



Since {2p + X)dww2 = F + {P{p) - P{p)), integrating by parts and using (pJj) . ([ST]) . ([3TT1) 
and ()3.66p . we find that 

h+l2<C (||P(p) - Pip)\\h + l|V^^|li2 + \\'^W2\\l2) 

+C\\P{p) - P{p)\\l4u\\l4^F\\l^ 
< \\\P^^^W2\\l2 + C {\\Vw2\\l2 + llVnIlia + Co) , (3.68) 

and using (j3.7p . (|3.1ip and (j3.66p . one easily gets 

h<C (^j p\w2\^d3^ (^J p\u\^dx^ llVtusllie 

< CC',°^'\\p'/'w2h-2 {Wp'^'mh^ + \\P{P) - PmLe) 

< C2Co'"/'||p^/2«;2|li2 + CC'J^. (3.69) 
Thus if Co is chosen to be such that 

Co<e5,2 = min{e5,i,(4C2)-'/^«}, 

then we obtain after putting (13.68j) . (13.69j) into (I3.67P and integrating it over (0,(t(T)) that 

sup ||Vt(;2||i2 + / \\ p^^'^ 71)21112 dt 

0<t<o-(T) Jo 

< CCo'/^ + C / (||Vw;2||i2 + \\Vu\\l2) dt < CCl'\ (3.70) 

Jo 

where we have also used (13. lip . (I3.60p and the following simple fact (due to (13. lip ) : 
|/o| < ^^\\Vw2\\l2+C\\P{p) - PiMh < ^\\Vw2\\l2+CCo. 
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We now estimate HVtysH^^a. To do so, using (|2.1|) . (|3.7|) . (|3.15|) and the standard H"^- 
regularity of elliptic system, we first deduce from (j3.58p that 



I|V«;3||l6 + ||V2u;3|| 



L2 

< C {Wpw^Wl^ + WVdV^WL^ + \\VvV''d\\L2) 

< c (wp'^'mh^ + c^'/'wvMl^ + INM'i'll'^M^'lNMlL^) ■ (3.71) 

Then, multiplying (|3.58|) by w^t in and integrating by parts, we find 
{filViusl"^ + {p + X){divw3f) dx + j p\m\'^dx 
^ ^ [Vd Vt> - \vd : Vvh] : Vw^dx 



1 d 
2di 



dt J \ 

- j ^Vd Vvt - ^Vd : Vvth^ : Vw^dx 

- j (vdt QVv- ^Vdt : V-uIs j : Vw^dx + J pw^ ■ {u ■ Vw3)dx 
d ^ 

i=l 

It readily follows from (j2.ip and p.lSp that 

|Jo| < CWVdUAWvU.W^w^WL^ < ^\\Vw3\\h+CCf'/'\\V\\\l,. (3.73) 
By (j3.15p . (j3.7ip and Cauchy-Schwarz inequality, we have 
Jl<C\\Vdh4Vvt\\L2\\Vw3\\Le 



< l\\p'^'n^3\\h+C{\\yvt\\h + W^Mh + llV^Hli^llV^dlli.) , (3.74) 
and similarly, after integrating by parts and using Cauchy-Schwarz inequality, we have 



J2<Cy {\dt\\V^v\\Vws\ + \dt\\Vv\\V^W3\)dx 

<C(MdlL2||V2^;||^3||Vt«3||L« + NdlL2||V7;|U^||VV||L0 

KCWdth^^M'/'WvM'/'Wv^wsh^ 
<c\\dth4^M]!'\NM]!' 



< \\\p'^'m\\h + C (llV^Hli. + llV^^lli.llV^dllt. + \\VMU\dt\\t.) ■ (3.75) 



Finally, similar to that in p.72p . we infer from (j3.7p . (|3.15p and p.7ip that 

c I I / p\u\^dx I 



\ 1/2 / /■ \ 1/3 

J3<C{ I p\w3\^dx] { p\u\^dx] \\Vw4l^ 
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< C,c'fyp'/'^,\\h + C {\\VMI2 + llV^^lli^llV^dlli.) . (3.76) 
Thus if Co is chosen to be such that 

Co<e5 = min{e5,2,(4C3)-3/^o}, 
then by (IXT2I1 and (f3:28]l we deduce after putting into ^^TTIb that 

sup ||Vi.3|li2 + / \\p"^W^\\l2dt 
0<t<a{T) Jo 

<C sup ||V2^;||2, +C / {\\Vvt\\l2 + \\VMh)dt 

0<t<cr(T) Jo 



0<t<cr{T) 

ra{T) 



+c sup \\vMh / {WdtWl. + dt 

0<t<a{T) Jo 

<C\\V^vo\\l2, (3.77) 
where we have used (|1.3p . p.l2p and (I3.15P to get that 



r wdtwi^dt <c r {\\v^dU2 + \\uvdu2 + \\vd\\i,fdt 

Jo Jo 



< 



T 



c / (I|v2(i||i2 + \\Vu\\L2\\vd\\L^ + wvdWL.W'^^dWL^ydt 







<C r {\\VH\\\2 + \\Vu\\\2)dt<CCl'\ 
Jo 



Similarly, multiplying (|3.72|) by t, integrating it over (0, o'(T)), and taking (|3.73|) - (j3.76|) into 
account, we deduce from (j3.29p and (j3.6ip that 

sup {twvw^wl^) + / tWp^l^w-iWl^dt 

0<t<a(T) Jo 



<C sup {t\\V^v\\l2)+C t{\,Vvt\\l2 + \\V\\\l2)dt 

0<t<a{T) Jo 

+C sup {t\\VMh) m\\- + \\^^d\\\-)di + C \Wn,Adt 

0<t<a{T) Jo Jo 

<C||V^;o|li2. (3.78) 
By the same token as that in the proof of (|3.65p , we conclude from (|3.77p and (|3.78p that 

sup ft'-f^\\Vws\\l2) + t'-^\\p'/'ws\\l2dt<C\\Vvo\\l,. (3.79) 

0<t<cr{T) ^ ''Jo 

Now, choosing wio = uq and vq = do so that wi + 11)2 + 103 = u and v = d, we immediately 
obtain (1333)1 from (l335]l . (ICTIl and ^l9\i . 
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To prove ()3.54p . choosing m = 2 — /? in (I3.46p . (I3.50p and adding them together, we infer 
from (I32D, (|3ll]) . IK2T\i and (l333]l that if Cq < £5, then 

sup t'-'(\\p'/'u\\h + \m\\h 

0<t<a{T) ^ 



+ / t'-^{\\V^l, + \\du\\l. + \\V^dt\\l,)dt 
Jo 

<c wvuWl^dt + c t^-^iivuii^^iiv^dii 

Jo Jo 



+c / t^-^ (^P^'^uWh + livdtiii^j dt + c f'^wvuWl.dt 

<C + C t^'I^WVulll^dt, (3.80) 







since it follows from (pm|) and (IX^ that for /3 G (1/2, 1], 

t'-'' (liv^lli. + livdiii.) [Wp'/'uWI. + \\vd,\\l, + wv'dwl,) dt 

< sup (t'~^{\\Vu\\l + \\VM\l,] 

0<t<a(T) ^ 

X i^-^ (iip'/'t^iii^ + iivdtiii. + iiv=^d|ii.) dt 

< c. 

Using Lemmas [m US (fXTTT) . (fXT2]) . (f3:2T1) and (l333]l . we have 

/ t^-^\\Vu\\l^dt<C t^-I^WVuWLiWVuWledt 
Jo Jo 
,a(r) 



<c£ t'^^WVuh^ (1 + Wp'/'uWh + \\V'd\\%'\\V^d\\%') dt 
<C + C (^ii-/3||vn|||,)'^' (i'-'^llp'/'^lli^) 







<C + C sup f(t^-^||Vn||i.)'^Vt^-/^||pV2^||i.)'^') r i^-'^llpV^^lli.dt 

0<t<(7(T) ^ ^ ' J Jo 

+c sup (t^-^liv^dlli.) M / ||v<.dt / t^-'^liv^dlli.dt 

0<t<<T{T) ^ ' \J0 ) \Jo ) 

<C + C7 sup [t^-^\\p^''^iifu^^'^ , (3.81) 

0<t<o-(T) ^ ^ 

provided Co < £5. As a result, putting ()3.8ip into ()3.80p and using Young's inequality, we 
immediately obtain p.54p . The proof of Lemma 13.61 is therefore completed. □ 
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We are now in a position of estimating A^^aiT)). 

Lemma 3.6 Let {p, u, d) he a smooth solution of f fj. j|) - |777P on M'^ x (0, T] satisfying ([57 
Then there exists a constant eg > 0, depending only on /j,, \, Aj'j, p, p, /3, Mi and M2, such that 

A,{aiT)) ^ sup / p\u\'dx < (3.82) 

0<t<a{T) J ^ 

provided Cq < Eq. 

Proof. Multiplying (jl.2p by 3\u\u and integrating it by parts over x (0, a{T)), we have from 
Cauchy-Schwarz inequality that 



sup / p\ufdx< I pQ\uQfdx + C I I \P{p) - P{p)\\u\\Vu\dxdt 
0<t<a{T) J J Jo J 

l-a{T) f- i-u(T) r 

+C j J \u\\Vu\^dxdt + C J J \Vd\'^\u\dxdt 
= J po\uofdx + "^Ri. (3.83) 



First, it is easily seen from ([2l]), dSZD and (l3TT]l that 

Ri<C \\P{p) - P{p)\\L3\\u\\Le\\Vu\\L2dt <C \\Vu\\\2dt < CCq. 

Jo Jo 



(3.84) 



Keeping in mind that /3 G (1/2, 1] and 5o G (0, 1/9], using Lemmas EH 1221 ([SZD, dXTB . (l3J2]l . 
p.lSp and Holder inequality, we deduce 



R2<C \\u\\L-o\\Vu\\l2dt <C \\Vu\\%^\\Vu\\2edt 

Jo Jo 

<C r WVuff {\\pu\\l2 + \\VdV^d\\L2 + \\Pip) - Pip)\\L<^f^ dt 
Jo 



< C r \\yu\\%^ {\\p^^^u\\l^ + \\^d\\LB \\V^d\\L2 + Cl^^^ ^ ^ dt 











<CCf''^ + C sup (a(3-W4||Vn||2,) 



5/4 



£ ^(2/3-3)/2^^ j l£ ^(3-2;3)/4 (^||pl/2^||2 ^ + WV^l^) dt\ 

< CCl^"'^, (3.85) 



24 



and similarly, 



Rz<C \\Vu\\^l^\\Vu\\Yi\\Vd\\l4V^d\\l2dt 
Jo 



U \\Vu\\l2dt\ U WV^dWl^dt 

+C sup f(a(3-W4||Vn||2,)'^V^{3-2W4||v2^||2 
0<t<a{T) \^ ^ ^ 

X ^^""^ ^3(2/3-3)78 (^^{3-2/3)/4 (^||^1/2^||2^ + ||v3d||2^ 



This, combined with I^SM^ . <K85\i and ([3:83]) gives 

sup [ p\u\^dx < [ po\uo\^dx + CCQ^°^^ 

0<t<a{T) J J 



< CC;^(2/3-l)/(4/3) ^ (^(^35o/2 < ^^(-35o/2^ ^g gg^ 



where we have used the interpolation and Sobolev embedding inequalities to get that 

\ 1/3 / f \ (2/3-l)/{4/3) / ^ X (3-2/3)/{12/3) 

Po\uofdx] <Ci po\uo\^dx\ i / \uof^^^^^^Ux] 

< CCl^'^-'^/^'^^WuoWf;^^ < CCf-'^'^'^^ < CCf'^'\ (3.87) 
since it follows from ()3.8p that 3(5o/2 < (2/3 — l)/(4/3). So, if Co is chosen to be such that 



C7o<e6=min{e5,(2Ci)"2A'}, 



then (|3.82p immediately follows from (j3.86p . □ 
We can now close the estimates of Ai{T) and A2{T). 

Lemma 3.7 Let (/>, u, d) he a smooth solution of \l.l\ )- [r7i^ on ]R3 x (0, T] satisfying \3. 7\ ). 
Then there exists a constant ej > 0, depending only on p, X, A,'y, p, p, 13, Mi and M2, such that 

Ai{T) + A2{T)<cl'\ (3.88) 

provided Cq < £7. 

Proof. To prove (fSrHSD . we first utihze Lemmas EH [221 ([XTT]) and (IXT^ to get that 

||F|li4 + ||6.||i4 <C{\\F\\l2\\VF\\1, + 

<c (\\Vu\\l.+cI'^) (\\p''^^h + C'^l^^d\\] 



lL2 ) ' 
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from which it follows that 



j a' [\\F\\% + dt 
Jo 

<C I a'^\\Vu\\L2\\p^/^u\\l2dt + C&o° [ a^\\Vu\\L2\\V^d\\l,dt 
Jo Jo 

^C&J'' r aV/^uWl^dt + CCl^'^'' ra^WVmidt 
Jo Jo 

4 

= ^Ri. (3.89) 



i=l 



In view of the definitions of Ai{T) (i = 1, . . . , 5) in (j3.ip - (|3.5p . we deduce from (j3.7p . (|3.53p 
and Holder inequality that 

Rl<C ^"^^^(6/3-7)78 (^^(3-2/3)/4||^^||2^y/2 (^^2||^1/2^||2^^ (^^1-/3 ||^l/2^||2^y/2 



lL2 



a(T) 



1/2 / \ 1/2 



{6/3-7)/4^^ 



+CAy\T)Al^\T) r a\\p'/'u\\l,dt 

Ja(T) 

< + CCo < CC^^^^°^'\ (3.90) 

Similarly, due to (|3.7p and (|3.2ip . we have 

4 „j' 



Y,R,< CC',^ I {a\\Vu\\l,f' {A\y'd\\l.f" {a'-^WV'dWl,) dt 

1=2 

A\p"'u\\hf' {a\\p'"u\\l.)dt 

+CCl"^'^ j\a'\\V'd\\l,f'\a\\V'd\\l,)dt 
Jo 

< C (c^''Ay\T)Al^\T) + cI^^Ai{T)aI^\t) + cl^^^^" Ai{T)AI^\t) 



< CCl^^^^°. (3.91) 
Thus, substituting (^Mij, PM]) into gives 

r (ll^llt^ + Mh) dt < C&i^''^'^ (3.92) 
Jo 

In terms of the effective viscous flux F, we can rewrite p.35p as 

= {P{p) - P{p))t + u ■ V{P{p) - P{p)) + -iP{p)diYU + 7(P(p) - P{p))dWu 
= (Pip) - P{p))t + u ■ V{P{p) - P{p)) + ^P{p)divu 
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+j;^iP(p) - nP))F + l^^iPip) - Pip)?, (3.93) 



which, multiphed by 2>a'^{P{p) — P{p))'^ and integrated by parts over M.^ x (0, T), yields 

< Ca'WPip) - Pip)\\h +C [ [cr'WPip) - Pmh + ^'llVnlli. + a^\\F\\l,) dt 

Jo 

< Cd^^^'^'\ (3.94) 

where we have used (I3.11|) . (I3.92|) and Cauchy-Schwarz inequahty. 

As a result of p.92p and ()3.94p . we deduce from the standard L^-estimate that 

a'^\\Vu\\\idt <C [ (||divu||^4 + ||a;||^4) dt 
Jo 

<C f a' (llFlli^ + \\P{p) - P{p)\\\. + Mt.) dt 

JO 

< CC'i^^'^'\ (3.95) 
Thus, it fohows from (pmi) . (pl^ . (pIMD and (1?:^ that 

^i(r) + A2{T) < CCo/^+^^° + C7Cj^+^°^/^ < C7iC7j^+^°)/^ < C^^^ 
provided Cq is chosen to be such that 



C7o<e7 = min{e6,C-'/'°}. 



The proof of Lemma 13.71 is therefore completed. □ 
Finally, we need to prove the uniform upper bound of the density. 

Lemma 3.8 Let [p, u, d) be a smooth solution of lll.l\ )- [T7T\ ) on x (0, T] satisfying |g. 7\ ). 
Then there exists a constant eg > 0, depending only on p, X, Aj'j, p, p, /3, Mi and M2, such that 

7 _ 

sup IIpIIloo < -p (3.96) 
o<t<T 4 

provided Cq < es- 

Proof. Let Dt — dt + u ■ \7 denote the material derivative operator. Then, in terms of the 
effective viscous flux F in p.3p . we can rewrite (II. ip as 

Dtp = g{p) + b'{p) 

where 

2p + X 2p + X Jq 
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Thus, to apply Lemma 12.31 we now need to estimate b{t). To do this, we first use ()2.2p and 
2.7p to deduce that for any < ti < t2 < (^(T), 



C / ||A-Miv(pn) + A-Mivdiv(M(d))||Loodt 
Jo 









lo 










= c 






'0 






< c 






'0 






< c 






'0 




+ i?2 



/.cr(T) 

Jo 



< C 







On one hand, using (|3.53p . (|3.54p and (|3.88p . we find 

Ri<C f^^\^^-'^/V/'n\\]!,' {t'-^\\Vn\\l,f" dt 

^(/3-2)/3||^l/2^||2/3^^\ Ij ^2-/3||^^||2^^^j 

£ ^-[(2-/3)(2/3-5o)+5o] (^^2-/3||^l/2^||2^y/3-^o (tUp^/'t^Hi^) dtj 

/ \ 3(l-<5o)/4 / (5.) X 35o/4 

<CU i-[2{2-/3)+3(/3-l)5o]/(3(l-5o))^^ j f j t||//2^||2 ^rft J 

since it follows from ()3.8p that 

2(2 - /3) + 3(/3 - l),5o ^ 9/3 -4/3^ + 1 
3(1 - 5o) ~ 7/3 + 1 

On the other hand, it follows from (j3.15p and (j3.2ip that 

R2 < CCf^''' j^^^^\-2il-P)/Z (^^l-/3||^3^||2^^2/3^^ 



> -1. 



1/3 / X 2/3 



Thus, putting the estimates of i?i,i?2 into p.97p gives 

\b{t2)-h{t^)\<C&^^'\ 
So, for t £ [0, (t(T)] one can choose Nq,Ni and in Lemma 12.31 as follows: 

iVo = CC5°/', 7Vi = 0, C = P- 
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Since it holds that 



we thus conclude from (12.91) that 



sup WpWl^ < max{p, p} + Nq < p + CCj'/^ < -p, (3.98) 

0<t<a{T) 2 



provided Cq is chosen to be such that 

Co < £8,1 - min <{ £j 



/ p \3/'5o 
V2cJ 

On the other hand, it follows from Lemmas EH [221 (jSTf]) . (f3lT]l and ^M\i that 

for any t G [cj(T),T]. Hence, using (|3.22p and (IS.SSh . we deduce that for any ti,t2 £ [a{T),T], 

|6(*2) - b{h)\ < C r \\F\\L^dt < 7r4-T(*2 -h) + C r \\Ff/ldt 
Jti Zp + A J^iT) 

+ f (a^Vu\\l.+a''l^\\V^dUA\'^'^d\\l,)dt 

2/i + A J^iT) V / 



Thus, for t G [(t(T),T] one can choose Nq,Ni and in Lemma 12.31 as follows: 

A 

2/i + A' 



Afo = CCo'/', iVi = ,^A_, r = /5 + L 



Noting that 



we can thus apply Lemma [2^ to get 



sup IIpIIl- < max<' + ij +iVo < — + CCo^^ < ^p, (3.99) 



o-(T)<t<T 

provided Cq is chosen to be such that 



Co < £8 = min |e8,i, | 



Therefore, the combination of (j3.98p with p.99p ends the proof of Lemma 13.81 □ 
Now, by virtue of Lemmas I3.1H3.8I we can complete the proof of Proposition 13.11 
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Proof of Proposition \3.1[ By Lemmas 13. 2( 13. 6113. 8| to prove Proposition I3.H it remains to 
estimate the term Ai{a{T)). In fact, using (j3.2ip . (|3.53p and (|3.88p . we have 

Ai{a{T)) < sup Icr ||Vn||^2 I sup (cr||Vn||^2j 



0<t<(7(T) ^ ^ 0<t<(7(T) 

/ i_fl,,„2 rii2 \ (2/3+l)/{4/3) 1,^2 7112 \ (2/3-l)/(4/3) 

+ sup [a^ f^WV^dWli) sup (crllV^dlll^a) 

0<t<cr(T) ^ ^ 0<t<o-(T) 

/ ..(T) \ (2/3+l)/(4/3) / (^) X (2,3-l)/(4/3) 

/ (^T) \ (2/3+l)/{4/3) / X (2/3-l)/{4/3) 

+ (y a^-'^llV^dlliadtj U a\\V^d\\l,dt\ 



provided Co is chosen to be such that 

Co<e^min{e8,(2C)-^2/3/(2/3-i)|_ 

Therefore, the proof of Proposition 13.11 is completed. □ 



4 Time-dependent higher-order estimates 

In this section, we prove the global estimates on the spatial-time derivatives of (p, u, d) 
which are needed to guarantee the existence of classical solutions. For this purpose, we assume 
that the conditions of Theorem 11.11 hold and that (p, u, d) is a smooth solution of (ll.ip - (|1.7p on 
X (0,T] satisfying Proposition 13.11 Moreover, from now on we will always assume that the 
initial energy Co satisfies (j3.10p . For simplicity, throughout this section we denote by C or Cj 
(i = 1, 2, 3, . . .) the various positive constants which may depend on 

I^A,Al,P,P,^,Mi,M2,g and T, 



where g £ L'^ is the function in the compatibility condition ()1.9p and T > is the time. 
First, one easily deduces from Lemmas l3.3H3.5l that 

Lemma 4.1 Assume that the conditions of Theorem hold. Then for any given T > 0, there 
exists a constant C(T) > such that 

sup {\\V^d\\l, + llVulliO 

0<t<T 

+ £ [Wp'^^iiWh + \\Vdt\\l2 + llV^dll^^) dt < C{T) (4.1) 

and 



sup [\\p'^'u\\i2 + WV'dWi^ + \\Vdt\\i2 

0<t<T ^ 
l-T 



+ I (l|Vu||i2 + \\du\\l2 + llV^dtllia) dt < C{T). 
Jo 



(4.2) 
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nil refit 



Proof. Indeed, choosing /3 = 1 in (I3.2ip and (I3.53p . one immediately obtains (14. 1|) from ()3.2ip . 
(I333]l and ^Mh . 

In order to prove (14.2p . we first notice tliat (I3.52p . togetlier with (14. ip . implies 

liv^diii^ <c(i + i|vdtiii2). 

Taking this into account, choosing m = in (|3.46p . (|3.50p and integrating them over (0, T), we 
have by ()4.ip and the compatibility condition ()1.9p that 

sup (llp^/^tillia + llV^dlll^ + \\Vdt\\l2 
0<t<T ^ 

+ [ (l|V^i||i2 + WdttWh + llV^dfllia) dt 

Jo 

<C + C [ \\Vu\\l4dt<C + C [ \\Vu\\l2\\Vi 
Jo Jo 

<C + cj^ {\\p^/^u\\l2 + W^d\\l2)dt, (4.3) 

where we have also used ()2.ip . (j2.6p and (j3.15p . Thus, in view of ()4.ip and (|4.3p . an application 
of Gronwall's inequality immediately leads to (j4.2p . □ 

Next, we use the Beale-Kato-Majda-type inequality (see Lemma [2. 4p to derive the estimates 
on the spatial gradient of the density as well as the velocity, which are very important for the 
proof of the higher-order estimates of the solutions. The proofs are mainly motivated by the 
ones in [TTlfTO]. 

Lemma 4.2 Assume that the conditions of Theorem \l.l\ hold. Then for any given T > 0, there 
exists a constant C{T) > such that 

sup (||Vp||i2ni6 + ||Vn||j,i)+ / \\Vuf/ldt < C{T). (4.4) 

Proof. For 2 < p < 6, it is easily derived from (jl.ip that 

^llVplliP < C\\Vu\\l^\\Vp\\lp + C\\V^u\\lp 

<C{1 + \\Vu\\loo) WVpWlp + C{1 + WpuWlp) , (4.5) 

where we have used (|4.ip . (j4.2p and the L^'-estimate of elliptic system to infer from ()1.2p that 

\\V^u\\lp < C {WpuWlv + ||V/>||lp + \\Vdy^d\\Lp) 

<C{\\pu\\lp + \\Vp\\lv + 1)- (4.6) 

We are now in a position of estimating ||Vu||loo. By (|2.ip . (j4.ip and (j4.6p . we deduce from 
(f2TT]l with g = 6 that 

||Vu||lcx) < C + C (||divti||j;,oo + ||V X u||loo) In (e + ||/0?i||L6 + ||Vp||^6) 
< C + C {\\(ii\u\\L°° + ||V X In (e + ||Vtt||/,2) 

+C(||div'u||Loc + ||V X n||Loo)ln(e + llVpll^e). (4.7) 
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So if we set 

llVplUe, ^{t) = 1 + (||divn||Loo + ||V x u\\loo) In (e + ||Vn||i2) + ||Vn||i2, 
then it follows from (|4.5p with p = Q and (|4.7p that 

< C7^'(f)$(t)ln$(t), 
due to the fact that ln<^(T) > 1. This particularly implies 

^ln$(t) < C^{t)\n^{t). (4.8) 
By virtue of Lemmas [HI [121 gH) and g^]), we find (w = V x u) 

\XjOC -\- W^WfjOo) dt 



< 



< 



^^/\t)dt <C + C [ {WViiWl^ + ||divn||^ 
Jo 

C + C I {\\P{p)-P{~p)\\l^ + \\F\\l^ + M\l^)dt 
Jo 

C + C [ {\\F\\l2 + \\uj\\l2 + \\VF\\l, + dt 
Jo 

<C + C I {\\Vu\\l2 + II Vdllioo \\V^d\\l2) dt < c, 
Jo 

and consequently, it follows from ()4.8p and Gronwall's inequality that 



(4.9) 



sup ||Vp(t)|U6 < sup ^{t)<C. (4.10) 

0<t<T 0<t<T 



As a result of ([421), and KWh . we obtain 

/" WVuf^^dt <C + C [ ^^/^{t)dt < C. (4.11) 
Jo Jo 

Using (14. ip and (14. lip , we also infer from (14. Sp with p = 2 and Gronwall's inequality that 

sup ||Vp(t)||i2 <C, 

0<t<T 

which, together with (|4.2p and (|4.6p . yields ||Vii||//i < C. The proof of Lemma 14.21 is therefore 
completed. □ 

Basing on Lemmas 14.11 and 14. 2^ we easily obtain 
Lemma 4.3 Assume that the conditions of Theorem \l.l\ hold. Then for any given T > 0, 

sup Wp'/^utWl^ + r \Wut\\l2dt < C{T), (4.12) 

0<t<T Jo 

sup {\\Vp\\Hi + \\VP\\m)+ [ \\VMHidt<CiT), (4.13) 

0<t<T Jo 

sup (llptll^i + llPtll^i) + r {WpttWh + WPttWh) dt < C{T). (4.14) 

0<t<T Jo 
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Proof. First, due to ut = u — u ■ Vn, one easily gets (I4.12p from Lemmas l2.lt I4.11 and 14.21 
Next we prove (|4.13p . To do this, noting that P{p) = Ap^ satisfies 

Pt + u • VP + 7Pdivu = 0, (4.15) 

from which and (jl.ip we have by some direct computations that 

Idiv^iii. + iiv^piii.) 

< C\\Vu\\loo (||VVlli2 + \\V^P\\1,)+C\\V^u\\l2 (||VV||l2 + ||V2p||^2) 
+C||V\||^6 (||Vp||i3 + ||VP||i3) (IIvVIIl^ + liv^piii^) 

< C (1 + \\Vu\\l^) {\\V^p\\l2 + WVPWh) + CllV^nll^^i, (4.16) 

where we have used ([XT]) and Using ()^7T|) . (021), (jHH), (|4.12p and Lemma we deduce 

from (|1.2p and the standard L'^-estimate of elhptic system that 

llV^nllj^i < C [WpiiWu^ + ||VP||j^i + \\SJ(N'^d\\H^) 

<C{l + \\Vut\\L^ + \\V^P\\L2). (4.17) 

Thus, combining (|4.17p with (j4.16p and using Gronwall inequahty, we immediately arrive at 
(jil^ since it follows from (03]) and (01^ that 

W^uf/l + \\'^ut\\l2 e Li(o,r). 

Finally, it is easily seen from (jl.ip and (j4.15p that 

WptWu^ + \\Pt\\m < C\\u\\l^ (||Vp||j^i + ||VP||j^i) + C\\Vu\\hi 
+C\\Vu\\L^{\\Vphs + \\VP\\L^) 
<C{1 + ||Vu|||i + llVpll^i + ||VP||^i) < C, (4.18) 

where we have used (j4.4p . (j4.13p and Lemma |2. II Moreover, noticing that (j4.15p implies 

Ptt + Uf'^P + u- VPt + 7Ptdivn + 7Pdivnt = 0, 

and hence, using (jl3D, (|4.12p . (|4.13p and (|4.18p . one obtains 

/ ||Ptt|li2dt<C [ {\\ut\\L4"^P\\L^i + \\u\\wl,^\\Pt\\Hl + \\Vut\\L2f dt 

Jo Jo 

<C + C [ {\\Vu\\l, + WVutWl^) dt < C. (4.19) 
Jo 

Analogously, one also has ||ptt||L2 G L^(0,T). So, combining this with (j4.18p . ()4.19p completes 
the proof of gUD. □ 



Due to the weaker compatibility condition (II. 9p . the methods used in |T7] to derive the 
higher-order estimates on the solutions {p,u,d), which guarantee the solutions obtained are 
indeed a classical one away from the initial time, cannot be applied any more. To overcome this 
difficulty, we need some careful initial-layer analysis. 
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(4.20) 



Lemma 4.4 Let a = min{l,t}. Then for any given T > 0, 

sup fT(||VM||^2 + llVutll^a + llV^dtll^a + llV^dll^a) + / \\V'^d\\jj2dt 
0<t<T Jo 

+ a (\\p''^uu\\l2 + \\Vut\\]j, + WVduWl^) dt < C{T). 

Proof. Differentiating (jl.2p and (jl.Sp witli respect to t gives 

putt - fJ.^ut - (/i + A)Vdivut = -ptut - {pu ■ Vu)t - VPt - divM{d)t (4.21) 

and 

dtt - Adt = {\Vd\^d)t - (n • Vd)t. (4.22) 

Tlius, multiplying (j4.2ip and (j4.22p by uu and —Adu respectively, and integrating the resulting 
equations over M^, we obtain after adding them together that 

Ijt j (/^l^"il^ + (^ + ^)(divut)^ + |Adt|2)(ix + j {p\utt\'^ + \Vdtt\'^)dx 

= - j ip^u, + . V. + pu, .Vu + pu. Vu, + VP, + divM(.),) . uud. 

+ j {ufVd + u-Vdt- I Vdpdi - 2d{Vd) ■ {Vdt)) ■ Adudx 

= -J^ -pt\ut\^ + {ptu-Vu + VPt + d\wM[d)t) -ut dx 

+ / (pun . Vu + p,u, . Vn + p,u • Vu, + d.vM(d),) • u^dx 

+1 / puM^dx - J Pudi^u^dx - J ipu, .Vu + pu. Vn.) • uudx 

+ j [{ut - Vd + u-Vdt - \Vd\'^dt - 2d{Vd) ■ {Vdt)) ■ Adu] dx 



i=l 



We now estimate each term on the right-hand side of (I4.23p . Using (jl.ip and integrating 
by parts, we first deduce from Lemmas 12. H I4.1H4.3I and Cauchy-Schwarz inequality that 

Io = - J [pu - {ut ■ Vut) + ptu ■ Vn • Ut - Ptdivut - M{d)t : Vnt] dx 
< C\\u\\L^\\p^^'^ut\\i,2\\Vut\\L2 + C||pt||2,2||n||L'^||Vn||i3||nf 11^6 

+ C\\Pt\\L2\\Vut\\L2 + C\\Vd\\L^\\Vdt\\LA\'^Ut\\L^ 

Since it holds that M{d)tt < C{\Vdtt\\Vd\ + jVfijp), we have by Lemmas EH [31H33] and 
integration by parts that 

h < C\\ptt\\LA\AL^\\^ALA\ut\\Le + C||/Ot||L2||^t|li6||Vn||i6 
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+C||/5t||i6||n||i6||Vnt||^2ht||i6 + C||Vd«||i2||Vd|U^||Vnt||i2 
+C\\Vdt\\]!,^\\Vdt\\%^\\Vuth2 

< \\\^du\\h + c {WpuWh + livt^illi. + Wv'dtWl.) . 

Due to (II. ip . one has pu = —div{pu)t, and hence, it follows from from (j4.12p . (|4.14p and 
integration by parts that 

I. + /3 = / ip^u . Vn, . + pu, . Vn, • d. - J PudWu,d. 

< C {\\pt\\L4u\\L^\\ut\\Le + \\pUt\\L4ut\\L6 + WPtth^) ||Vnt||i2 

< C [WVutWh + llp'/'n,||i/'||n,||^/2||Vn,||i2 + \\Pu\\h) 
<C{l + \\Vut\\l, + \\Pu\\l2). 

and similarly, by (|4.4p one gets 

h <C{\\ut\\L4^u\\L3 + ML^WVuth^) Wp^^^uWl^ 
<l\\p'^'uu\\h+C\\Vut\\l,. 

Finally, integrating by parts and using Lemmas 12.11 and I4.1H4.31 we can estimate the last 
term on the right-hand side of (I4.23P as follows: 

/5 < {\Vut\\Vd\ + IntllV^dl + \Vu\\Vdt\ + \u\\V^dt\) \Vdu\dx 

+c j {\yd\\y^d\\dt\ + ivdpivdti + \y^d\\ydt\ + \yd\\y^dt\) \ydtt\dx 

<C||Vd||i«.||Vnt||i2||Vd«||i2+C7||nt||i6||V2d||^3||V<i«||i2 
+C7||Vn||i3||Vdt||L6||Vdtt||i2 +C7||n||L-||V2dt||i2||Vdtt||i2 

+ C||Vd|U^||V2d||i3||dt||,.6||V<i«||,.2+C7||Vd||ioo||Vdt||i2||Vd«||i2 

+c\\v^d\\L4Vdt\\LA\^dtt\\L^ + c||vd||L-||v2dt||i2||vdtt||i2 

< \\\^du\\l. + C (1 + ||Vnt||i2 + WV^dtWl,) . 

By virtue of the estimates of /j (i = 0, 1, . . . , 5), we deduce after multiplying (j4.23p by cr{t), 
integrating it over [0, T] and using Gronwall's inequality that 

sup a {\\Vut\\l2 + llV^d^ll^,) + r a (Wp^/^ttWh + llVd^tHi^) dt < C, (4.24) 

0<t<T Jo ^ ' 

where we have also used (fO|) . (|4.12p and (|4.14p . This, together with (|4.13p and (|4.17p . gives 

sup (a||Vtx||^2) < sup a(l + ||Vnt||22 + ||VP|||i) <a (4.25) 
o<t<r o<t<T 

On the other hand, by Lemmas 14 . Ifl4 . 3 1 we obtain by applying the standard L^-estimate to 
(OB that 

||V2ut||i2 < C\\p^ut + (^ + A)Vdivnt||i2 



35 



< C\\ {put + pu-Vu + VP + divM{d))^ \\l2 

< cWp'/'uuWl^ + cWpth'AWth'^ + cWptWML^WVuhe 

+C\\ut\\Le\\Vu\\L3+C\\u\\L^\\Vut\\L^+C\\VPt\\L2 
+ C\\Vd\\Loo\\V^dt\\L2 + C\\Vdt\\L4^^d\\L3 

<c(l + \\Vut\\L2 + Wp'/^uuWl^ + WV^dth^) , (4.26) 



which, together with (|4.24p . yields 

f a\\V'^ut\\l2dt < C. (4.27) 
Jo 

In a similar manner as the derivation of ()4.25p . we also infer from l\1.3\i that 

< C {WV^dtU^ + \\V\u . Vd)\\L2 + \\V\\Vd\''d)\\L2) 
<C{1 + WV^dtU^ + \\Vu\\h^ \\Vd\\H2 + livdii^.) 
<C + C\\V^dt\\L2, 

from which, (|4.24p . and (02]), it follows that 

sup {a\\V''dfH2)+ r \\V^d\\]j2dt<C. (4.28) 

0<t<T JO 

Therefore, combining (gJl]), (jOSj) . (jOT]) and ()i:28]) finishes the proof of ()T20]) . □ 

To prove the Holder continuity of the first-order derivatives of density and pressure, we need 
the following lemma which is concerned with the VF^''?-estimate (g G (3,6)) on the gradients of 
density and pressure. 

Lemma 4.5 For fixed q £ (3, 6), it holds for any T > that 

sup (IIVpIIh/1. + IIVPIIh^,,) + r {\\Vutrj^> + IIV^^II^,,) dt < CiT), (4.29) 

0<t<T JO 

where 

1 <Po < e (1,2). (4.30) 

oq — o 

Proof. Operating to both sides of (jl.ip . (j4.15p . and multiplying them by g|V^p|'^~^V^/0, 
q\V'^P{p)\'^^^'S/^P{p), respectively, we obtain after integrating by parts over and using (|4.4p . 
(j4.13p and Lemma [27T] that for any q E (3,6), 

I (llvVllL + llv^PiiL) 
< c\\vuu^ (llvVllL + \\y'P\\h) + ciiv^nii^i,, (llvVlli;' + liv^pii^;^^ 
+c\\vMl^ (lIvpllL-llvVlli;' + \\vpu^\\v'p\\r/ 

<Cil + \\Vu\\H2){l + \\V'p\\% + \\V'P\\%) 

+c\\vMw^^^ f llvVlli;' + llv^Piii;^ . (4.31) 
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Applying the standard L^-estimate to the elhptic system (|1.2p yields that 

\\V^u\\wi,, < C\\V^u\\li + C\\V{put + pu-Vu + VP + divM{d))\\Li 

< C\\Vu\\h2 +C{\\Vp\\L4ut\\L^ + llVntlliO 

+C {\\u\\l^\\V'^u\\li + ||Vn||Loo ||Vn||L9 + ||V/?||l9 ||Vu||l=o) 

+C7||V2p||i, + C (||Vfi||Loo||V='(i||L, + \\V^d\\Loo\\V^d\\L',) 

< C (1 + \\Vut\\L2 + WVutU', + \\Vu\\h2 + \\V^d\\H2 + ||v2p||i,) , (4.32) 
where we have also used Lemmas 12.11 and I4.1H4.31 Putting (|4.32p into (j4.3ip gives 

^ (IIVVIIL + IIV'PIIL) <C{1 + \\Vuth2 + WVuth. + \\Vu\\j,2 + \\V''d\\H2) 

X (l + ||vVlli, + ||W||i,). (4.33) 
For po being the one in (|4.30p . it foUows from (|4.12p and (|4.20p that 

r {\\Vut\\L2 + \\Vu\\h2 + \\V^d\\H2y' dt 
Jo 

<C + C snv ia\\Vu\\]j2+(j\\V^d\\]j2Y°''^ a'^'^l'^dt < C, (4.34) 

0<t<T Jo 



and moreover, using Lemma |2. 11 Holder inequality and (j4.20p . we find that 

Jo 



^ n f -Po/2 / ||V7 ||2 \Po(6-(j)/(4g) / ,,2 

Jo 

(\ Po(6-g)/(4g) 

IIV7 Il2 I / -pn/2 /' IIV7 Il2 \Po(3g-6)/{4(5r) 

0<t<T / Jo 



dt 

f fT \ (4g-po{3g-6))/(4g) . ^ \ Po(3g-6)/(4g) 

<cf / a-2fo''/(4'?-Po(39-6))^^\ M ^iivntll^idtj 

< C, (4.35) 

since < 2pog'/(4g - ^0(30^ - 6)) < 1 and < po(3g - 6)/(4g') < 1. 

By virtue of (|4.34p and (j4.35p . we deduce from (j4.33p and Gronwall inequality that 

sup (||V/9||h/i,, + ||VP||iyi,,) < C, Vge(3,6), 

0<t<T 

which, together with (|i32D . and (iOSjl . yields the desired estimate of (g^S]). □ 

Finally, we need the following initial-layer analysis which particularly implies (uf, V^n) are 
Holder continuous away from the initial time t = 0. 

Lemma 4.6 For any given T > 0, it holds that 

sup a (\\p'^/^utt\\L2 + ||VV||i2 + ||V\||h/i,,) + / a^\\Vutt\\l2dt < C{T). (4.36) 

0<i<T ^ ^ Jo 
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Proof. Differentiating (I4.21|) with respect to t gives 

puttt + pu ■ ^utt - pAutt - ip + X)Vdivutt 

= 2div{pu)utt + div{pu)tut - 2{pu)t ■ Vut - (pttu + 2ptUt) ■ Vu 
-puu ■ Vu - VPtt - dwM{d)u, 

wliich, multiplied by u^t in and integrated by parts over M^, yields 
/ + J {P'Nuttl'^ + {p + X){divuttf) dx 

= -^jpu. Vuu ■ uud. - j\pu), . (V(.. . + 2Vu, . d. 
-jipuu + 2p,u,).Vu.uud.-jpuu-Vu.uud. 

+ I Pttdivuttdx + / M{d)tt : Vuttdx = ^ (4.37) 

i=i 

The right-hand side of (j4.37p will be estimated term by term as follows, using Lemma |2.H 
I4.1H4.5I and Cauchy-Schwarz inequality as well. 

Ji < \\u\\l^\\p^^^uu\\l4'^uu\\l^ <6\\Vuu\\l2+C{6)\\p^/\tt\\l2, 

J2 < C{\\put\\L3 + ||^i||L°°||pt||L3) (l|Vnt||i2||ntt||^6 + ||ui||i6||Viiit||i2) 



<6\\Vuu\\l2+C{5){l + \\Vut\\l), 

<5\\Vnu\\l2+C{6){\\pu\\h + \\Vut\\l.), 
J4 < C\\p'/\tt\\L4^u\\Ls\\uu\\L^ < 6\\Vuu\\l2 + C{6)\\p^/^uu\\l2, 
J5 < C\\Pu\\L4^uu\y < 6\\Vuu\\l2+C{6)\\Ptt\\l2 

and 

Je <c(\\Vd\\L^\\Vdu\\L2 + \\^dt\\]!,^\\Vdt\\%-) WVutth^ 

< 6\\Vuu\\l2 + C{6) {\\Vdu\\l2 + llV^dtllia) , 

since direct computations give 

\M{d)u\<C{\Vd\\Vdu\ + \Vdt\^) . 

Thus, putting the estimates of Ji, . . . , Jg into (j4.37p and multiplying the resulting inequality 
by cr^ , we have by choosing 6 > small enough and using Lemmas 14.11 - 14.41 that 



sup (a^\\p^^^utt\\l2) + / a^\\Vuu\\l2dt 

<t<T ^ ^ JO 

<C + C r a\\p^/^uu\\l2dt + C r {\\Ptt\\l2 + \\ptt\\l2 + a\\Vdu\\l2)dt 
Jo Jo 
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+C sup a^/^\\Vut\\L^ + \\V^dt\\L2) r {\\Vut\\l2 + \\V^dt\\l2)dt 

0<t<T Jo 

< C. (4.38) 
As a result of ()4.38p . ()4.26p and Lemma 14.41 we also see that 

fT||VV||L2 <C(j(l + \\Vut\\L2 + Wp^/^uuWv^ + WV^dtU^) < C, (4.39) 
and thus, it follows from (j4.32p . (I4.39p . Lemmas 14.41 and [451 that for any q G (3,6), 

(t\\V^u\\w^,, <Ca[l + ||Vnt||i2 + ||Vnt||L, + ||Vn||j^2 + \\V^d\\H2 + \\^^P\\li) 
<Ca{l + WVutWui + ||Vu||h2 + ||V^d||j^2 + \\V'^P\\li) 
< C. (4.40) 

Therefore, combining (j4.38p - (|4.40p immediately leads to ()4.36p . □ 



5 Proof of Theorem 11.11 

With all the a priori estimates at hand, we are now ready to prove our main results. To 
this end, we first need the following local existence theorem of classical solutions of (|l.ip -( fL7|) 
with large initial data. 

Proposition 5.1 Assume that the initial data {po,uo,do) satisfy the conditions il.8\) . il.9\) of 
Theorem \l.l\ Then there exist a positive time Tq > and a unique classical solution {p,u,d) of 
(QHP-ir^j on X (0,To], satisfying p>0, \d\ = 1, and for any r G (0,To), 

({p-p,P{p)-P{p)) G C{[0,To];H^nW^''i), 
{p - p, P{p) - P{p)) G L-(0, To; H TV^.a) 
u G C([0, T]-D^r\ D^) n L°°(r, Tq- n D^''i), 
utGL^{T,To-D^nD^)nH\T,To-D^), 
Vd G C([0, To];H^) n L-(t, Tq; H^), 
[ dt G C([0, To];H^) n L°°(t, Tq; H^). 

Proof. As that in |20], we can use the Galerkin's approximation method to construct the 
approximate solutions to the momentum equation, then use this approximate and the 
equations of conservation mass and angular momentum to get p"^ , . The existence of a smooth 
approximate solution (/o™, n'", d*") follows from the fixed point theorem, similar to that on the 
compressible Navier-Stokes equations (see, for example, |20p36j). Now, in order to prove the 
convergence for the approximate solutions and to obtain a smooth solution of ()l.ip - ()1.7p . it is 
essential to derive some uniform a priori estimate for (p^ , u"^ , d"^) . 

It has been shown in [20| Theorem 2.1 or (3.21), (j3.26p ] that there exists a small Tq > 0, 
independent of m and the lower bound of density, such that 



(5.1) 



sup {WV^uTWh 

0<t<To 
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+ ihniU, + llV^lli, + ||V^d-||i, + llV^dnii.) < Cexp (Cll^lli.) , (5.2) 

where C > may depend on Tq, but is independent of the size of domain. 

In view of the bounds in (15.2p . we can proceed to derive more higher-order estimates on the 
solutions (/o™", n*", d*") in the same way as those carried out in Lemmas l4.3H4.6l To summarize 
up, for any r G [OjTq] the approximate solutions {p^ ,(F^) satisfy 



sup iWVp^Ww^. + ||vp-||p^i,o + / ° {WVuTWZ + liv^«"^||^ij dt 

<t<To Jo 

+ sup {\\^u^\\l, + \\v\n\m^w^. + \N<\\m)+ Tin 

T<t<Tn Jt 



+ sup [\\y''d'^\\]j, + \\V^dT\\l2)+ r\\VdTt\\l2dt<C (5.3) 

T<t<To Jt 

for some positive constant C which may depend on Tq, but is independent of m and the lower 
bound of density. 

With the help of the local estimates ()5.2p and ()5.3p , we easily deduce after taking a subse- 
quence [p^^ , vT^^ , d"^^ ) and passing to limit as j — t- oo that (p^^ , , d^^ ) would converges to a 
solution {p, u, d) of (|l.ip - (|1.7p on M'^ x (0, Tq) satisfying (|5.ip due to the lower semi-continuity. 
The uniqueness of strong/classical solutions can be proved in the same manner as that in |21] . 
This finishes the proof of Proposition 15.11 □ 

Proof of Theorem By Proposition 15.11 there exists a small time Tq > such that the 
Cauchy problem (ll.ip - (ll.7p has a unique classical solution {p,u,d) on x (0,To]. We shall 
make use of the a priori estimates, Proposition 13.1! and Lemmas I4.lfl4.51 to extend the local 
classical solution {p, u, d) to all time. 

First, in view of the definitions in (I3.ip - (l3.5p . it is easily seen from (I3.20p and (13.87P that 

Ai(0)+^2(0) = 0, ^3(0)<Co^°, ^4(0) + ^5(0) <Co^» and p^ < p, 

due to Co < e. Thus, there exists a Ti G (0,To] such that (13.71) hold for T = Ti. 
To be continued, we set 

T* = sup{r I dST]) holds} . (5.4) 

Then T* > Ti > 0. Hence, for any < r < T < T*, it follows from ([Ll]), gai), (gH]), ([TOD . 
(fOnjl . (I06D and Lemma ED that 

j j \dt{p\ut\^)\dxdt + j j \dt{p\vSJu\^)\dxdt 

< C y j {\pt\\ut\'^ + p\ut\\utt\) dxdt 

+C j j (|pt||n|2|Vn|^ + p|n||nt||V'up + p|'up|Vu||V'Ui|) dxdt 

< C y j [p\(ii\u\\ut\'^ + \u\\V p\\ut\'^ + p\ut\\utt\) dxdt 

+C j j (|divti||up|Vnp + |Vp||u|^|Vnp + |np|Vn||Vut| +^1^11^411X7^^) (ixdt 
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<C J (^||Vu||Loo||/3^/2ut||^2 + ||w||L6||V/9||i2||ut||^6 + \\p^^'^Ut\\L2\\p^^'^Utt\\L2^ dt 

+C j (||Vn||f^2 + ||Vu||f^2||Vp||i2 + ||Vn||f^2||Vnt||i2 + \\p^^^ut\\L2\\Vu\\%2') dt 
<C{r,T), 

which, together with (|4.4|) and (|4.12|) . yields 

p'/^ut, p'/^u-Vu€C{[t,T];L^). 

and consequently, 

p'/'ueC{[T,T];L^). (5.5) 

Next, we claim that 

T* = oo. (5.6) 

Otherwise, T* < oo. Then by Proposition l3.1l (|3.9p holds for T = T*. So, it follows from Lemmas 
HTTOl and (|53D that {p,u,d){x,T*)) satisfies (1.11) and ([L9]), where g{x) = p^/'^u{x,T*). Thus, 
Proposition 15.11 implies that there exists some T** > T* such that (|3.7p holds for T = T**, 
which contradicts (j5.4p . Hence, (j5.6p holds. Proposition 15.11 and Lemmas I4.1H4.6] thus show 
that (p, n, d) is in fact a unique classical solution on x (0, T] for any < T < oo. 

In order to complete the proof of the existence part, it remains to prove that (p, u, d) is 
continuous in t, especially, to prove that 

{p-p,P{p)-P{p))(^C{[G,T\■D^r^D^^'i), <7G(3,6), (5.7) 

since by virtue of Lemmas I4.1H4.5I one easily deduces from (ll.l|) ~ (ll.3p that 

{{p-p, Pip) - Pip)) e C([0, T];H^n W^'"), 
ip -p, Pip)- Pip)) eCi[0,T];H^nW^''^- weakly), (5.8) 

^iu,Vd) e C([o,r];i?i nz)2). 

To prove (15. 7p . we denote by Dij = dfj with i,j = 1, 2, 3. Then it follows from (II. ip that 

dtDijp + div(nDjjp) = —divipDiju) — div(9i/) • dju + • diu) 

holds in D'(IR'^ x (0,T)). Now, let i;/(x) is the standard mollifying kernel with width u and set 
p'^ = p* jy. Then, we infer from the above equation that 

dtDijp" + diviuDijp") = -divipDiju) * jy - div(9ip • djU + djp ■ diu) * ju + Ru, (5.9) 

where R^, = dxviuDijp^) — diviuDijp) * jy satisfies (cf. [HU Lemma 2.3]) 

f \\RA\%r.L.dt <C f \\u\Q,^^\\D,,pr;\^^^dt < C (5.10) 
JO Jo 

due to (|4.29p . where po > 1 being the same one as in ()4.30p . 

Multiplying ()5.9p by q\Dij p'^\'^~'^ Dij p'^ and integrating by parts over M^, we obtain 



d ur. 

-rJDijp' 



dt 



\\Li 



-iq—1) J {Dijp'^l'^divudx — q J divipDiju) * ju\Dijp'^\'^ '^Dijp'^dx 
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-q j (div(5ip • djU + djp ■ diu) * jy) \Dijp''\'^ "^Dijp^dx 



which, combining with ()4.29p and (jS.lOp . yields 



sup ||VV||l9 + 
o<i<r 

<C + C 
< C. 



2 ^i^UQ 



dt 



This, together with Ascoh-Arzela theorem, gives 

||VV"(-,OIIl^ ^ l|VV(-,i)llL« in <^([0,T]) as ^ 0. 
In particular, we have 



Similarly, one also has 



|v^p(-,t)IU, GC([o,r]). 



|vV-,i)llL2 GC'([o,r]), 



(5.11) 



which, combined with ()5.1ip and (j5.8p 9. implies 

vVg C([0,r];L2nL'?). 



(5.12) 



In the same way, we can also prove that V2p(p) G C([0,r];L2nL''). This, together with (f5T2]l . 
finishes the proof of (|5.7p . and thus, the existence of a classical solution {p,u,d) of (jl.ip - (|1.7p 
has been proved. 

Next we prove the large-time behavior of {p,u,d). This can be done as the ones in |16pi7). 
however, for completeness we sketch the proof here. Multiplying (I3.35P by 4(P(p) — P{p))^ and 
integrating it over R^, we get that 

^^\\P{p) - P{p)\\h = J {\Pip) - P{p)\'diyu - 2,^P{p){P{p) - P{p)fdi^u) dx, 
which, together with (|3.94p and (|3.95p . shows 



±\\P{p) - P{p)\\l, 



/oo 
{\\Vu\\l, + \\P{p)-Pip)\\l^)dt<C. 



As a result, we have 



\\Pip) - Pip)\\L^ ^ as t^oo. 
This, together with (|3.1ip and the uniform upper bound of density, shows that 

lim \\p — p\\ip = 



(5.13) 



holds for any p as in (jl.l4p . 
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To study the large-time behavior of the velocity, we set 

M(t)^^||V^||i. + ^||div^||i.. 
Then, multiplying (ll.2p by ii in and integrating by parts over M^, we obtain 

- j (fiAu + + A)Vdivn) ■ iidx = j {{P{p) - P(p))divn + M{d) : Vu - dx. (5.14) 

Recalling the definition of " ' " , we deduce after integrating by parts that 



^^||Vn|||2 + /i j (^dkU^dkUidiU^ - ^\Vu\'^{d\vu)^ dx, 



and similarly, 

— /(// + A)Vdivu • iidx = — 
J 2 dij 



ifi + X)VdiYU ■ iidx = ^i±^i-||divw||i2 + OdiVnIlia), 



which, inserted into (j5.14p . yields 

\M'{t)\ < C (llVnIlis + \\p"^u\\l. + (||P(p) - + ||Vd||i3||Vd||i6) WViiU^ 

< C {\\Vu\\L2\\Vu\\\i + \\p^''^u\\\2 + ||Vu||i2 + ||V2(i|||2||Vn||i2^ . 

Thus, by (IHISHI) and (13:95]) we see that 

|M'(t)|2dt < cy (l|p^/^^lli2 + ||Vn||^4 + llVttll^a) (it 



< C sup 11/9^/2^ 



1/2,-, II 2 



) \\p^''^u\\hdt < C, 



,t>i 

from which and ()3.88p it follows that 

\\Vu(t)\\L2 ^ as t^oo. (5.15) 

As a result, we also have 

J p^/^\u\'^dx < (^J p\u\'^dx^ \\u\\l6 <C\\Vu\\l2 ^ as t ^ oo. (5.16) 

Moreover, using (jS.lip . (j3.15p and (|3.88p . we infer from (j2.6p that 

sup||Vn(t)||i6 < C, 
t>i 

which, together with (I5.15P and the interpolation inequality, leads to 

||Vn(t)||Lr ^ as t^oo, VrG[2,6). (5.17) 
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Finally, applying V to both sides of ()1.3p and taking the L^-inner product, by (I3.1ip and 
(j3.88p we deduce after integrating by parts over x (1, oo) that 



/oo 
{WVdtWl^ + \\V^d\\l, + \\Vd\\% + llVdlliallV^dllie) dt 

/oo 
(l|v^x|ii.||vd|iioo + hiiieiiv^diii^iiv^dii^e) 

/oo 
(llVfitllia + llV^dll^^ + \\V^d\\l, + llVnll^^) dt 



<C, 

which, together with p.lip . gives 

\\V^d{t)\\L2 ^ as t^oo. 
This, combining with (jS.lip and (j3.88p . yields that for any k £ (2,6), 

\\Vd{t)\\^^i,k ^ as t^oo. (5.18) 

Combining (|5.13p with (j5.15p - (|5.18p immediately proves (|1.13p . The proof of Theorem 11.11 is 
therefore completed. □ 
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